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What controls star formation in the central 500 pc of the Galaxy? 
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ABSTRACT 

The star formation rate (SFR) in the Central Molecular Zone (CMZ, i.e. the central 500 pc) of 
the Milky Way is lower by a factor of ^ 10 than expected for the substantial amount of dense 
gas it contains, which challenges current star formation theories. In this paper, we quantify 
which physical mechanisms could be causing this observation. On scales larger than the disc 
scale height, the low SFR is found to be consistent with episodic star formation due to secular 
instabilities or variations of the gas inflow along the Galactic bar. The CMZ is marginally 
Toomre-stable when including gas and stars, but highly Toomre-stable when only accounting 
for the gas, indicating that the condensation of self-gravitating clouds may be limited. On 
small scales, we find that the SFR in the CMZ is consistent with an elevated critical density 
for star formation due to the high turbulent pressure - potentially aided by weak magnetic 
effects and an underproduction of massive stars due to a bottom-heavy initial mass function. 
The existence of a universal density threshold for star formation is ruled out, as well as the 
CO . importance of the H1-H2 phase transition of hydrogen, the tidal field, the magnetic field, 

radiation pressure, and cosmic ray heating. We propose observational and numerical tests to 
distinguish between the remaining candidate star formation inhibitors, in which ALMA will 
play a key role. We conclude the paper by proposing a self-consistent cycle of star formation 
in the CMZ, in which the plausible star formation inhibitors are combined. Their ubiquity 
suggests that the perception of a lowered central SFR should be a common phenomenon in 
other galaxies. We discuss the implications for galactic star formation and supermassive black 
I hole growth, including a prediction that the recently reported bimodality of star formation in 

high-redshift galaxies may emanate from a difference in the gas inflow rates. 

Key words: Galaxy: centre — galaxies: evolution — galaxies: ISM — galaxies: starburst — 
galaxies: star formation — stars: formation 
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1 INTRODUCTION 

Star formation in galactic discs is often described with a power law 
relation between the star formation rate surface density Esfr and 
the gas surface density E (Kennicutt 1998b): 

Esfr. = AskE jv , (1) 

as was originally proposed for volume densities by Schmidt (1959). 
The typical range of power law indices is N = 1-2, whether 
E refers to dense or all gas, and across a range of spatial scales 
(Kennicutt 1998a,b; Bigiel et al. 2008; Liu et al. 2011; Lada et al. 
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2012; Kennicutt & Evans 2012). Star-forming discs in the local 
Universe follow a tight, 'Schmidt-Kennicutt' relation with Ask = 
2.5 x 10~ 4 and N = 1.4 ± 0.1 (Kennicutt 1998b), with Esfr in 
units of M© yr _1 kpc -2 and E in units of Mq pc~ 2 . When using 
only molecular gas, A mol = 8 x 10 -4 and TV" = 1.0 ± 0.2 (Bigiel 
et al. 2008, 2011). Since N > 1 indicates some effect of self- 
gravity, we refer to these relations as density-dependent star for- 
mation relations. Another commonly-used expression is the Silk- 
Elmegreen (Silk 1997; Elmegreen 1987, 1993, 1997) relation: 

Esfr = AseSJI, (2) 

with fi the angular velocity at the edge of the star-forming disc 
and Ase = 0.017 (Kennicutt 1998b) a proportionality constant, 
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Figure 1. Three-colour composite of the CMZ. with in red the HOPS NH3(1, 1) emission (Walsh et al. 201 1; Purcell et al. 2012) to indicate the gas with a 
density above a few times 10 3 cm -3 , in green the MSX 21.3/xm image (Egan et al. 1998; Price et al. 2001), and in blue the MSX 8.28/jm image. The MSX 
data shows PAH emission (mostly tracing cloud edges), young stellar objects, and evolved stars. The labels indicate several key objects and regions. 



adopting the same units as before and writing Q in units of Myr - . 
As for the case of N > 1 in equation (1), the dependence on the 
angular velocity implies that this star formation relation is density- 
dependent. 

These galactic-scale, global star formation relations have the 
advantage that they are easily evaluated observationally, but the de- 
pendence of surface densities on projection suggests that more fun- 
damental physics drive the observed scaling relations. Recent anal- 
yses of star formation in the solar neighbourhood have been used 
to propose a possibly universal, local volume density threshold for 
star formation riLada ~ 10 4 cm" 3 above which most 1 gas is con- 
verted into stars on a ~ 20 Myr time-scale (Gao & Solomon 2004; 
Heiderman et al. 2010; Lada, Lombardi & Alves 2010, although 
see Gutermuth et al. 201 1 and Burkert & Hartmann 2012 for an op- 
posing conclusion). It has been argued that this threshold also holds 
on galactic scales (Lada et al. 2012), in which case it would connect 
low-mass star forming regions to high-redshift starburst galaxies. 

It has been shown that Espr drops below the relations of 
equations (1) and (2) beyond a certain galactocentric radius (e.g. 
Martin & Kennicutt 2001; Bigiel et al. 2010) - the straightfor- 
ward detection of cutoff radii is well-suited for testing star for- 
mation thresholds. However, these are also the regions of galaxies 
where the minority of the dense gas resides. The central 500 pc of 
our Galaxy (the Central Molecular Zone or CMZ) contains a few 
per cent of the total molecular gas mass in the Galaxy and also a few 
per cent of the total star formation. At a first glance, this suggests 
that the star formation rate in the Milky Way is simply directly pro- 
portional to the reservoir of molecular gas available to form stars. 
However, the average density of molecular gas in the CMZ is two 
orders of magnitude higher than that in the disc. Specifically, the 
CMZ contains ~ 80% of the NHa(l, 1) integrated intensity in the 
Galaxy (Longmore et al. 2013a), reflecting an overwhelming abun- 
dance of dense gas (n > a few x 10 3 cm -3 ). These two facts are 
in direct contradiction with the proposed volumetric and surface 
density star formation relations that predict that a given mass of 
gas will form stars more rapidly if the density is higher (e.g. Kenni- 
cutt 1998b; Krumholz & McKee 2005; Padoan & Nordlund 2011; 
Krumholz, Dekel & McKee 2012). In Longmore et al. (2013a), we 
show that these relations do indeed over-predict the measured star 

1 Except for the mass loss due to protostellar outflows (e.g. Matzner & 
McKee 2000; Nakamura & Li 2007). 



formation rate by factors of 10 to 100. The low SFR in the CMZ is 
particularly striking because its gas surface density is similar to that 
observed in starburst galaxies, which seem to follow a sequence 
with a factor of ten higher SFRs than predicted by typical star for- 
mation relations (Daddi et al. 2010b). If the SFR depends on the 
density, something is required to slow the rate of star formation in 
the CMZ compared to the rest of the Milky Way and other galaxies. 

The SFR in the CMZ is also inconsistent with the Lada et al. 
(2012) star formation relation, despite the fact that it does not rely 
on the density. Because most gas in the CMZ is residing at densi- 
ties larger than the Lada threshold, a gas consumption time- scale 
of 20 Myr implies a SFR of 1 M@ yr -1 , which is 1-2 orders of 
magnitude higher than the measured SFR (Longmore et al. 2013a). 
In nearby disc galaxies, a simple proportionality of the SFR to the 
molecular mass is commonplace (Bigiel et al. 2008, 2011). De- 
spite the similarities between the CMZ and high-redshift galaxies 
(Kruijssen & Longmore 2013), which do form stars at or above 
the rate predicted by density-dependent star formation relations, the 
CMZ is consistent with an extrapolation of the molecular star for- 
mation relation with a constant H2 depletion time-scale (i.e. N = 1 
in equation 1, also see below). This is surprising - self-gravity im- 
plies that dynamical evolution proceeds faster when the density is 
higher (i.e. N > 1). If gravity is an important driver of star for- 
mation in the CMZ and galaxy discs, then a constant molecular gas 
depletion time-scale requires that some resistance to the gravita- 
tional collapse towards stars must increase at a rate comparable to 
self-gravity when the density goes up. Due to its extreme charac- 
teristics, the CMZ is the prime region to understand the underlying 
physics. 

In this paper, we take the point of view that starbursts and 
other regions which suggest TV > 1 (and thereby imply a density- 
dependence to the SFR) do not seem to apply to the CMZ in the 
Milky Way. To understand this difference better, we evaluate the 
global and local processes that affect the rate of star formation in 
the central region of the Milky Way. Thanks to the strongly con- 
trasting environments of central and outer galactic regions, galaxy 
centres provide a unique and independent way to study the univer- 
sality of star formation relations. We exploit this contrast to discuss 
the implications of the lack of star formation in the CMZ for exist- 
ing star formation relations with N > 1. The paper is concluded 
by sketching a possible picture of how local and global star forma- 
tion criteria connect, and we propose observational and numerical 
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Table 1. Properties of the considered regions and galaxies. 



Object ID 


^CO,20 
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Qgas 


2*,2 


Qtot 


EsFR 


^depl 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(?) 


(8) 


(9) 


(10) 


(ii) 


CMZ 230 pc-integrated 




1.2 


0.75-1.2 


20-50 


0.5 


(9.0-36) 


38 


(1.2-2.1) 


0.20 


(0.6) 


CMZ 1.3° cloud 




2.0 


0.88-1.5 


20-30 


0.5 


(6.4-16) 


20 


(2.1-4.0) 


0.13 


(1.5) 


CMZ 100-pc ring 




30 


1.6-3.3 


25-50 


0.1 


(0.95-4.0) 


29 


(0.76-2.7) 


3.0 


(1.0) 


Kennicutt (1998b) low E disc 


2.8 


0.05 


(0.04) 


(2.6) 


(0.5) 


1.5 






0.0024 


(2.1) 


Kennicutt (1998b) high E disc 


2.8 


0.20 


(0.07) 


(5.9) 


(0.6) 


1.5 






0.017 


(1.2) 


Daddi et al. (2010b) low E starburst 


0.4 


3.0 


(0.34) 


(6.1) 


(0.05) 


1.5 






3.9 


(0.08) 


Daddi et al. (2010b) high E starburst 


0.4 


10 2 


(0.80) 


(86) 


(0.3) 


1.5 






560 


(0.02) 


Daddi et al. (2010b) low E BzK 


1.8 


2.0 


(0.05) 


(55) 


(6) 


1.5 






0.28 


(0.7) 


Daddi et al. (2010b) high E BzK 


1.8 


10 


(0.07) 


(200) 


(10) 


1.5 






2.7 


(0.4) 



-^CO,20 = Xco/10 20 (K km s" 1 cm 2 )- 1 ;» 0.53a C o/M© (K km s" 1 pc 2 )" 1 is the CO-to-H 2 conversion factor, E 2 = E/10 2 M pc~ 2 is the 
gas surface density, k is the epicyclic frequency in units of Myr -1 , a is the velocity dispersion in units of km s _1 , = h/10 2 pc is the scale height, 
Qgas is the Toomre (1964) disc stability parameter when only including the self-gravity of the gas, E*,2 = S*/10 2 Mq pc -2 is the stellar surface density, 
Qtot includes the self-gravity of gas and stars, Esfr is the star formation rate density in units of Mq yr _1 kpc~ 2 , and t dcpl = EEg FR /Gyr is the gas 
depletion time. Values in parentheses are implied by the other numbers (see text). 



tests through which the different components of this picture can be 
addressed in more detail. 

Throughout the paper, we adopt a mean molecular weight of 
= 2.3, implying a mean particle mass of /im H = 3.9 x 10~ 24 g. 



2 OBSERVATIONAL CONSTRAINTS FROM THE 
CENTRAL MOLECULAR ZONE 

Longmore et al. (2013a) present the observational constraints for 
the suppression of star formation in the CMZ, which are sum- 
marised here. A division is made between global and local physics, 
where 'global' refers to spatial scales larger than the disc scale 
height (AR > h), on which the ISM properties are set by galactic 
structure, and 'local' refers to spatial scales smaller than the disc 
scale height (AR < h), on which the ISM properties are set by mi- 
crophysics. For instance, the formation of giant molecular clouds 
(GMCs) proceeds on global scales, whereas processes internal to 
the GMCs, such as stellar feedback, the turbulent cascade, or mag- 
netic fields then determine the fraction of the GMC mass that pro- 
ceeds to star formation. 



2.1 Global constraints 

Figure 1 shows a three-colour composite of the central 4° of the 
Milky Way, corresponding to a spatial scale of ~ 600 pc at the 
distance of the Galactic centre (we adopt 8.5 kpc, which is con- 
sistent with Reid et al. 2009). We combine the dense gas emission 
(NH3(1, 1), in red) with infrared imaging (green and blue) to high- 
light the gas close to (or above) the Lada, Lombardi & Alves (2010) 
density threshold, young stellar objects, and evolved stars. The leg- 
end indicates the objects in the CMZ that are most relevant to this 
paper. In particular, the 1.3° cloud and the 100-pc, twisted ring of 
dense gas clouds between Sgr C and Sgr B2 (Molinari et al. 201 1) 
will be scrutinised in detail. The GO. 253+0.016 cloud (the 'Brick') 
was recently identified as a possible progenitor of a young massive 
cluster (Longmore et al. 2012), and we will use it as a template for 
the massive and dense clouds that populate the CMZ. 

In Table 1 , we list the derived properties of the gas in the CMZ, 
as well as the typical characteristics of galaxies as defined by exist- 
ing, empirical star formation relations. The columns indicate (1) the 
object ID, (2) the CO-H2 conversion factor used to derive the gas 



mass, (3) the gas surface density, (4) the epicyclic frequency, (5) 
the velocity dispersion, (6) the scale height, (7) the Toomre (1964) 
Q stability parameter of the gas disc (see §3.1.1), (8) the stellar sur- 
face density, (9) the Toomre Q parameter of the entire disc, (10) the 
observed star formation rate surface density, and (1 1) the gas deple- 
tion time-scale. The first three rows list the properties of the CMZ, 
where we distinguish two components (the 100-pc ring described 
by Molinari et al. 2011 and the 1.3° cloud), as well as its proper- 
ties smeared out over a radius of 230 pc (corresponding to 1.55°, 
i.e. just beyond the 1.3° cloud). This division into sub-regions is 
made because the 100-pc ring and the 1.3° cloud contain most of 
the dense (and thereby total) gas mass within \l\ 1.5°. 

The total gas mass within 230 pc is taken from Ferriere, 
Gillard & Jean (2007), whereas the observed gas surface density 
in the plane of the sky of the 100-pc ring and the 1.3° cloud is 
derived from the HiGAL data (Molinari et al. 2010) following the 
analysis in Longmore et al. (2013a). Note that this estimate relies 
on dust emission and is insensitive to the conversion of CO inten- 
sity to H2 column density (the ratio of which is represented by the 
parameter Xco), and hence is not affected by any uncertainty in 
Xco (cf. Sandstrom et al. 2012, and references therein). However, 
the derived column density does depend on the assumed gas-to- 
dust ratio, which may be lower in the CMZ than in the Galac- 
tic disc due to the higher metallicity (see Longmore et al. 2013a 
for a discussion). At the high densities of the CMZ, almost all of 
the gas is in molecular form. The 3D geometry of the 100-pc ring 
was taken from Molinari et al. (2011, assuming a ring thickness of 
10 pc), whereas the face-on projected geometry of the 1.3° cloud 
as an ellipse with semi-major and minor axes of 85 pc and 50 pc 
is taken from Sawada et al. (2004). We adopt the rotation curve 
of Launhardt, Zylka & Mezger (2002, Figure 14) to calculate the 
epicyclic frequency. 2 The velocity dispersion and scale heights of 



2 Note that this gives circular velocities of 140-200 km s _1 for the three 
regions listed in Table 1 , whereas the peak line-of-sight velocity of the 1 00- 
pc ring as measured from the HOPS NH3(1,1) emission (Walsh et al. 201 1; 
Purcell et al. 2012) is 80 km s at the location of Sgr B2. Such a low 
line-of-sight velocity may be caused by the possible eccentric orbit of the 
ring (Molinari et al. 201 1), which is thought to lie under such an angle that 
the edge of the ring lies close to apocentre as seen from Earth. In this sce- 
nario, the measured line-of-sight velocity should be lower than the local 
circular velocity at the position of Sgr B2, a difference that would be ampli- 
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the dense gas in both the 100-pc ring and the 1.3° cloud are deter- 
mined from the HOPS NH 3 (1,1) emission (Walsh et al. 2011; Pur- 
cell et al. 2012). The value ranges in column 5 of Table 1 represent 
the maximum and minimum measured velocity dispersions across 
each region. We calculate <5 gas from the preceding columns as in 
equation (3) in §3.1.1 below. As can be seen in Table 1, the stel- 
lar component dominates the potential. We obtain the stellar sur- 
face density from the adopted rotation curve (Launhardt, Zylka & 
Mezger 2002), by deriving the spherically symmetric volume den- 
sity profile and only including the stellar mass in a slab of thickness 
of 2h, Therefore, the total stellar surface density of the bulge at the 
same galactocentric radii is higher than listed in Table 1. Column 9 
gives Toomre Qtot, which is corrected for the presence of stars as 
in equation (5) in §3.1.1 below, using a stellar velocity dispersion 
of a* ~ 100 km s" 1 (de Zeeuw 1993). Note that if we assume that 
the 1.3° cloud is roughly spherically symmetric, its virial parameter 
(Bertoldi & McKee 1992) a ~ 2 indicates it is roughly in equilib- 
rium, but only when including the stellar gravity. If the stars would 
be absent, the cloud would be highly unbound and hence stable 
against gravitational collapse. A similar effect of the stellar poten- 
tial is seen when comparing Toomre Q of columns 7 and 9 for the 
230 pc-integrated CMZ. By contrast, the Q parameter of the 100-pc 
ring is hardly affected by the presence of stars. The SFR densities 
Ssfr are derived using Table 2 of Longmore et al. (2013a), which 
lists the SFR for certain parts of the CMZ as derived from free-free 
emission (Lee, Murray & Rahman 2012). For details on the deriva- 
tion, we refer to the discussion in Longmore et al. (2013a). The gas 
depletion time follows from E and Esfr as £de P i = EEgp H . While 
the listed values reflect the depletion times of all gas, i.e. including 
both Hi and H2, it is important to note that most of the gas is in 
molecular form. Only for the disc sample from Kennicutt (1998b) 
this does not hold, where the molecular depletion times are lower 
than those for all gas listed in Table 1 . 

The fourth and fifth rows of Table 1 span the star formation 
relation of equations (1-2) for the disc galaxy sample used in Ken- 
nicutt (1998b). The bottom rows span the high-SFR 3 and low-SFR 
sequences of Daddi et al. (2010b). At a given surface density, the 
typical epicyclic frequencies of the Kennicutt (1998b) galaxies are 
obtained by combining equations (1) and (2), whereas for Daddi 
et al. (2010b) they follow from their equations (2) and (3). The ex- 
tragalactic velocity dispersions are obtained by assuming Q — 1.5 
and using the definition of Q (see §3.1.1). The characteristic disc 
scale heights are added to distinguish between the aforementioned 

fied further by possible projection effects. An alternative explanation is that 
the 100-pc ring extends to higher longitudes than stated in Molinari et al. 
(201 1), in which case the circular velocity would exceed the observed line- 
of-sight velocities at the presumed tangent points due to projection alone 
(Bally et al. in prep.). This picture is consistent with the proper motion of 
Sgr B2, which is ~ 90 km s _1 (Reid et al. 2009). Using the measured 
velocities instead of the Launhardt, Zylka & Mezger (2002) rotation curve 
gives a factor of ~ 1.6 lower epicyclic frequencies. Both extremes are used 
to calculate the possible range of ft listed in Table 1, and hence also con- 
tribute to the range of Q ga s and Qtot • 

3 This extra sequence of star-forming galaxies contains some galaxies from 
the nearby starburst sample of Kennicutt (1998b), but has an elevated SFR 
with respect to that paper because a different value of Xqo is assumed 
(Daddi et al. 2010b). We emphasise that inferring the H2 density using CO 
emission is an indirect method and hence may introduce substantial un- 
certainty. However, recent results from the KINGFISH survey of nearby 
galaxies suggest that Xco weakly decreases with star formation rate den- 
sity (e.g. Sandstrom et al. 2012), which is at least qualitatively consistent 
with Daddi et al. (2010b). 
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Figure 2. Observed star formation rate density as a function of the global 
star formation relations from Kennicutt (1998b). Left: using the Schmidt- 
Kennicutt relation from equation (1), Right: using the Silk-Elmegreen rela- 
tion from equation (2). The blue symbols span the sequence of nearby disc 
galaxies from Kennicutt (1998b), the red symbols span the star-forming 
galaxies from Daddi et al. (2010b), the open and closed, black symbols in- 
dicate the 1.3° cloud and the 100-pc ring, respectively (see Table 1), and 
the grey symbol denotes the spatially integrated CMZ. The solid lines in- 
dicate the 1:1 agreement, and the dotted lines in the left panel represent 
Esfr = j4 m olX! (bottom, using A mo i = 8 X 10 — 4 as in Bigiel et al. 
2008) and £ S fr = ^skS 2 (top). 

'global' and 'local' regimes, and are calculated by assuming an 
equilibrium disc including a factor-of-two increase of the disc self- 
gravity due to the presence of stars (cf. Elmegreen 1989; Martin & 
Kennicutt 2001). 

We compare the data from Table 1 to the global star formation 
relations of equations (1) and (2) in Figure 2. The spatially resolved 
elements of the CMZ (open and closed, black symbols) are forming 
stars at a rate that is a factor of 3-20 (i.e. typically an order of mag- 
nitude) below either relation. By contrast, the Schmidt-Kennicutt 
relation does describe the CMZ well when spatially smoothing it 
over a 230 pc radius (the closed, grey symbol), whereas smoothing 
hardly affects the agreement with the Silk-Elmegreen relation. It is 
not obvious whether this simply implies that the global star forma- 
tion relations break down at spatial scales smaller than ~ 500 pc. 
Smoothing is justified on galactic (J; 500 pc) scales, because the 
lifetime of substructure on these scales is typically of the order of 
(or shorter than) an orbital timescale, and accounting for substruc- 
ture therefore only introduces spurious stochasticity (e.g. Schruba 
et al. 2010). However, the nuclear rings that appear in numerical 
simulations of barred galaxies are persistent over many dynamical 
times (e.g. Piner, Stone & Teuben 1995; Kim, Seo & Kim 2012), 
and hence it seems physically incorrect to smear out the 100-pc 
ring and any other, possibly persistent structure in the CMZ to a 
much larger scale when comparing to global star formation rela- 
tions. These structures should be accounted for when describing the 
star formation (or lack thereof) in the CMZ. Therefore, the physi- 
cally appropriate representation of the CMZ in Figure 2 is given by 
the open and closed, black symbols. 

Also visible in the left-hand panel of Figure 2 is that the 
CMZ (which is nearly completely molecular) agrees well with the 
Bigiel et al. (2008) star formation relation between the SFR and 
the molecular gas mass (shown as the lower dotted line), which 
was derived for nearby disc galaxies with surface densities 3 < 
Eh 2 /Mq pc -2 < 50. The CMZ extends this range by roughly two 
orders of magnitude. Bigiel et al. (2008, 2011) show that the gas 
depletion time of the galaxies in their sample is td cp i = 1-2 Gyr, 
which is indeed consistent with time-scales listed for the CMZ in 
Table 1. In the framework of this relation, the CMZ is somehow the 
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norm while both the Daddi et al. (2010b) samples of BzK (Tacconi 
et al. 2010; Daddi et al. 2010a) and starburst/submillimeter galax- 
ies (Kennicutt 1998b; Bouche et al. 2007; Bothwell et al. 2010) are 
the exception. The Bigiel relation is the only global star formation 
relation that fits the CMZ - the region remains anomalous in the 
context of equations (1) and (2), or the Lada et al. (2012) relation 
(see §2.2). As will be shown in §3, there are several possible rea- 
sons why the CMZ is peculiar. 

2.2 Local constraints 

Turning to local physics, in Longmore et al. (2013a, Figure 4) we 
have shown that 70-90 per cent of the gas in the CMZ resides at 
column densities above the Lada column density threshold, and 
hence should be forming stars. Throughout this paper, we use the 
corresponding volume density threshold of nLada ~ 10 4 cm~ 3 
(Lada, Lombardi & Alves 2010). While this may apply for the low 
(~ 10 cm ) densities of GMCs in the solar neighbourhood, the 
mean volume density of the gas in the CMZ is no ~ 2 x 10 4 cm~ 3 
(Longmore et al. 2013a), i.e. more than two orders of magnitudes 
higher, and it hosts GMCs with typical densities of ~ 10 5 cm -3 . 
As mentioned in §1, the threshold for star formation and deple- 
tion time-scale from Lada et al. (2012) would imply a SFR that is 
1-2 orders of magnitude higher than the measured value (Long- 
more et al. 2013a). If we assume that the SFR is driven by self- 
gravity, we can use the observed SFR to derive the fraction of the 
gas mass residing above some volume density threshold for star 
formation. Given that the free-fall time above nLada is ts < 1 Myr 
and that the observed depletion time-scale of molecular hydrogen 
in the CMZ is ~ 1 Gyr (see the last column of Table 1), we con- 
clude that if a local density threshold for star formation n t h exists, 
only 0ttff/e CO reide P i Si 0.5 per cent of the gas mass in the CMZ 
should be above this threshold. 4 

The density probability distribution function (PDF) dp/dn of 
a turbulent interstellar medium (ISM) is often represented by a log- 
normal, of which the width and dispersion are set by the Mach 
number (e.g. Vazquez- Semadeni 1994; Padoan, Nordlund & Jones 
1997). However, the high-density tail can also (at least locally) be 
approximated by a power law, i.e. dp/dn oc n~^ (e.g. Klessen 
2000; Kritsuk, Norman & Wagner 2011; Elmegreen 2011; Hill 
et al. 2012). This can be used to estimate the required value of n t h 
in the CMZ. As mentioned above, the fraction of gas that is used 
to form stars is a factor of /sfr,„ ~ 0.005 lower than the fraction 
above the Lada threshold. This reduction implies that for the as- 
sumption of a power law density PDF n th = /gp^~^7iLada- This 
relation is shown in Figure 3, which gives 4 x 10 8 > n th /cm~ 3 > 
10 5 for exponents 2.5 < /3 < 4. The density PDF of the gas in 
the CMZ is unknown, but the typical high-density slope due to 
self-gravity in the numerical simulations of Kritsuk, Norman & 
Wagner (2011) is f) = 2.5-2.75, which suggests nth = 10 7 - 
4 x 10 8 cm -3 . Note that because a power-law tail implies some 
effect of self-gravity, this functional form should fail to describe 
the detailed shape of the PDF in the low density-regime where self- 
gravity is not important. This needs to be kept in mind when choos- 
ing a roughly representative value of /3. 

4 This estimate assumes that the maximum star formation efficiency in 
protostellar Cores IS £core 

= 0.5 and the star formation timescale is 
4> t tg = 1.91% as in Kramholz & McKee (2005, Eqs. 19-21). However, 
the precise numbers are not important to the order-of-magnitude argument 
that is made here. 
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Figure 3. Volume density threshold for star formation implied by the lack 
of star formation in the Central Molecular Zone of the Milky Way as a 
function of the assumed power law slope of the density PDF at densities 
n > "Lada (solid line). The dashed line indicates the required threshold 
density for a log-normal PDF (Padoan, Nordlund & Jones 1997) with Mach 
number M = 70 and mean density no = 2 X 10 4 cm~ 3 , appropriate for 
the CMZ (see text). The dotted line indicates the same, but accounts for the 
effect on the density PDF of a magnetic field with strength B ~ lOOuG at 
a temperature of T = 65 K. 

For a log-normal PDF with a Mach number of M ~ 70 
(cf. Table 1, with a temperature of T = 65 K as in Ao et al. 
2013 - also see Morris et al. 1983) and a mean density of no = 
2 x 10 4 cm -3 , the required threshold density is even higher at 
n t h ~2x 10 9 cm~ 3 . However, this does not account for the in- 
fluence of the high magnetic field strength (B ~ 100uG, Crocker 
et al. 2010) near the Galactic Centre. For a thermal-to-magnetic 
pressure ratio 2cl/v\ < 1, with c s the sound speed and va the 
Alfven velocity, the dispersion of the log-normal density PDF is 
suppressed. Assuming T — 65 K, we find 2c$ /v\ — 0.31 for 
the CMZ. If we modify the density PDF accordingly (cf. Padoan 
& Nordlund 2011; Molina et al. 2012, adopting B oc n 1/2 as in 
Padoan & Nordlund 1999), the threshold density required by the 
observed SFR becomes n t h ~ 5 x 10 8 cm -3 . We conclude that 
the CMZ firmly rules out a density threshold at n = 10 4 cm -3 , and 
adopt a lower limit of nth = 10 7 cm -3 . This is still exceptionally 
high in comparison to the threshold density of nearby disc galax- 
ies, and is only known to be reached in dense, rapidly star-forming 
galaxies at high redshift (e.g. Swinbank et al. 201 1). 



3 POSSIBLE MECHANISMS FOR INHIBITING STAR 
FORMATION 

In this section, we summarize and quantify which physical mech- 
anisms may limit the SFR in the central regions of galaxies with 
respect to the SFRs predicted by density-dependent star formation 
relations. The implications of the candidate inhibitors for existing 
star formation relations are discussed in §4, where we also pro- 
pose ways of distinguishing their relative importance observation- 
ally. This section is separated into global and local star formation 
inhibitors. We discuss in §5 how these can be connected. 

3.1 Globally inhibited star formation 

We first discuss the possible mechanisms that may limit star forma- 
tion in the CMZ on spatial scales larger than the disc scale height 
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(AR > h). A lack of star formation may be caused by a stability of 
the gas disc against gravitational collapse, or by global dynamical 
processes that limit the duration of star formation - either stochas- 
tically or due to long-term dynamics. 



3.1.1 Disc stability 

A recurring question regarding galactic-scale star formation rela- 
tions has been in which part of the parameter space they apply - 
is there a threshold density below which star formation is negli- 
gible? And if so, by which factor is the SFR reduced? The ex- 
istence of such a threshold is suggested by the sharply truncated 
Hn discs in galaxies, which indicate that beyond a certain galac- 
tocentric radius Esfr falls off more rapidly than suggested by 
equations (1-2) (Kennicutt 1989; Martin & Kennicutt 2001; Bigiel 
et al. 2010). Note that a radial truncation is absent for the molecular 
(N = 1) star formation relation (Schruba et al. 2011). The possible 
physics behind surface density thresholds are extensively discussed 
by Leroy et al. (2008), but we briefly summarize them here. 

The efficiency of galactic star formation may be related to the 
global gravitational (in)stability of star-forming discs - if the kine- 
matics of a gas disc are such that it can withstand global collapse, 
star formation is suppressed (Toomre 1964; Quirk 1972; Fall & Ef- 
stathiou 1980; Kennicutt 1989; Martin & Kennicutt 2001). There 
are indications that the threshold density for gravitational instabil- 
ity not only applies globally, but also locally for the spiral arms and 
the inter-arm regions of M51 (Kennicutt et al. 2007). The density 
threshold for gravitational instability in galaxy discs is based on the 
Toomre (1964) Q parameter: 



Qgas — 



ttGE' 



(3) 



with a the one-dimensional gas velocity dispersion and k the 
epicyclic frequency, which is a measure for the Coriolis force in 
a condensing gas cloud: 



f-V / dlnV\ 



din VV' 2 



(4) 



with V the circular velocity, R the galactocentric radius, and fl = 
V/R the angular velocity. Note that for a flat rotation curve we 
have k — 1.41fl, whereas for solid-body rotation k = 2Q. In gas 
discs with Q < 1, the self-gravity of contracting clouds is sufficient 
to overcome the Coriolis force and undergo gravitational collapse, 
whereas Q > 1 indicates stability by kinetic support. Galaxy discs 
typically self-regulate to Q ~ 1 (e.g. Martin & Kennicutt 2001; 
Hopkins, Quataert & Murray 2012), with an observed range of 0.5 
to 6 for galaxies as a whole (e.g. Kennicutt 1989; Martin & Kenni- 
cutt 2001), and an even larger variation within galaxies (see below 
and Martin & Kennicutt 2001). If a substantial mass fraction of the 
disc is constituted by stars, the disc can be unstable even though 
Q gas > 1. In that case we write (cf. Wang & Silk 1994; Martin & 
Kennicutt 2001): 



Qtot = Q 1 + 



E a J 



(5) 



with cr* the stellar velocity dispersion. 

If star formation is driven by global disc instability, then the 
inversion of equations (3) and (5) leads to a critical surface density 
for star formation, modulo a proportionality constant uq: 



where ip = 1 if the contribution of stars to the gravitational poten- 
tial is neglected (in the solar neighbourhood tp ~ 1.4, see Martin & 
Kennicutt 2001). Kennicutt (1989) and Martin & Kennicutt (2001) 
empirically determined for nearby disc galaxies that ctQ = 0.5- 
0.85 with a best value of ctQ ~ 0.7, indicating that the criti- 
cal Toomre Q parameter for star formation is Qcrit = 1.2-2 or 
Qcrit ~ 1.4. The threshold density depends on o and k, both of 
which increase towards the Galactic centre (e.g. Morris & Serabyn 
1996; Oka et al. 2001; Longmore et al. 2013a; Shetty et al. 2012), 
and hence it is possible that gravitational stability inhibits star for- 
mation in the CMZ. 

A different, but closely related form of global disc stability to 
self-gravity is due to rotational shear, which may compete with self- 
gravity and prevent the collapse of the disc to form stars (Goldreich 
& Lynden-Bell 1965; Elmegreen 1987, 1991; Hunter, Elmegreen & 
Baker 1998). Rather than the epicyclic frequency of equation (6), 
this threshold depends on the shear time, i.e. the time available for 
gas instabilities to arise during the shear-driven density growth of 
spirals (Elmegreen 1987, 1993, 1997). This shear condition may be 
more relevant than the Toomre condition if the angular momentum 
of a growing gas perturbation is not conserved, as might be the case 
in the presence of magnetic fields or viscosity. The shear time-scale 
is the inverse of the Oort (1927) constant A: 



A ort = -0.5R^-. 

dR 

The critical density for self-gravity to overcome shear is then 

O"^4oort 



(7) 



(8) 



Sloomrc — CtQ 



nGtp irGQcvittp ' 



(6) 



Eoort = a a- , . 

nGtp 

with a.A ~ 2.5 (Hunter, Elmegreen & Baker 1998). Like k, A Qolt 
decreases with galactocentric radius, and hence shear motion may 
be responsible for the suppression of star formation in the CMZ. 

Before continuing, we should caution that the Toomre or shear 
stability of a gas disc can only act as a 'soft' threshold. There are 
several ways in which such a threshold may be violated, and as a 
result it represents a soft separation between star-forming and qui- 
escent gas discs. For instance, small-scale structure (e.g. spirals or 
bars), turbulent dissipation, magnetic stripping of angular momen- 
tum, or a soft equation of state could all lead to local star formation 
in a disc that is globally stable to star formation. 

The overdensities E/E{ ToomrCi0 ort} w i m respect to the sur- 
face density thresholds for gravitational instability and overcoming 
rotational shear are shown in Figure 4 as a function of galactocen- 
tric radius for a simple model of the Milky Way. For R < 0.3 kpc, 
we only include the three data points from Table 1, at galactocen- 
tric radii of R = {0.08, 0.19, 0.23} kpc for the 100-pc ring, 1.3° 
cloud, and the 230-pc integrated CMZ, respectively. The resulting 
overdensities are shown with and without the stellar contribution 
to the gravitational potential. For R 3 kpc, we use the Wolfire 
et al. (2003) model for the ISM surface density profile, including 
a factor of 1.4 to account for the presence of helium, and adopt 
a = 7 km s _1 (Heiles & Troland 2003). At these galactocentric 
radii, the epicyclic frequency k and the Oort constant A are calcu- 
lated using the rotation curve from Johnston, Spergel & Hernquist 
(1995). Our conclusions are unaffected when using other Milky 
Way models. Finally, we correct for the presence of spiral arms 
by dividing £{ Toomrc>0 ort} by a factor of two (cf. Balbus 1988; 
Krumholz & McKee 2005). 

Figure 4 shows that the Milky way disc is unstable to star for- 
mation in a ring covering 4 ^ i?/kpc ^ 17. Especially the outer 
edge of the star-forming disc catches the eye. The difference be- 
tween the Toomre and Oort thresholds across the disc is generally 
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Figure 4. Ratio of the gas surface density in the Milky Way £ to the critical 
surface density for star formation £{Toomrc,Oort} as a function of galac- 
tocentric radius. The solid and dashed lines refer to the critical densities of 
equations (6) and (8), respectively. The red lines include the stellar gravita- 
tional potential, whereas the blue lines exclude the contribution from stars. 
The range of Toomre Q implied by the ratio S/SToomro is indicated on 
the right-hand side, and the error bars indicate the range of uncertainty from 
Table 1 . Star formation in the Galactic disc is suppressed in the grey-shaded 
area, i.e. for R 4 kpc and R 20 kpc. 

too small to indicate with certainty which mechanism dominates. 
Only in the 100-pc ring we tentatively find Eioomrc > Soort, 
which indicates a decreasing importance of shear as the circular 
velocity decreases for i? < 0.1 kpc, but the difference is less 
than la. A more substantial change is brought about by includ- 
ing the presence of stars in the calculation of Q (red lines in 
Figure 4). When the stellar gravity is ignored, the CMZ outside 
of the 100-pc ring is highly stable to gravitational collapse with 
E/S{ ToomrCi oort} ~ 0.1, suggesting that star formation could po- 
tentially be suppressed. However, when the stars are included, the 
CMZ appears marginally stable - the range of Q to t measured in the 
CMZ is in fact very similar to that observed in normal disc galax- 
ies (compare Table 1 to Martin & Kennicutt 2001), which suggests 
a similar degree of self-regulation. The paucity of star formation 
in the CMZ is therefore not due to the 'morphological quench- 
ing' (Martig et al. 2009) of star formation that is considered to en- 
able the long-term presence of quiescent gas reservoirs in galactic 
spheroids if Qtot > 1. 

The high Q gas may still slow down the condensation of self- 
gravitating gas clouds and their decoupling from the stellar back- 
ground potential. The time-scale for clouds to become gravitation- 
ally unstable is i grav ~ Qgas/^ (e.g. Jogee, Scoville & Kenney 
2005). If the SFR is limited by slow cloud condensation, this there- 
fore implies a decrease of the SFR by 1/Q gas . This simple mod- 
ification of the Silk-Elmegreen relation is consistent with the ob- 
served SFR for the 230 pc-integrated CMZ, but it does not explain 
the other components of the CMZ. Note in particular that the 100- 
pc ring is always marginally Toomre-stable, both in terms of Q gas 
and Qtot, which is also consistent with its clumpy, beads-on-a- 
string morphology (Longmore et al. 2013b). We return to this point 
in §5.3. 



(i.e. a few Myr) that the current reservoir of dense gas may not be 
related to the observed star formation tracers, which originate from 
gas that was present at least one dynamical time-scale ago. The 
CMZ agrees with the Schmidt-Kennicutt relation when averaged 
over a size scale of R ~ 230 pc, which corresponds to an orbital 
period of 14 Myr (see Table 1). It cannot be ruled out that these are 
the size and time-scales on which the variability of star formation is 
sampled well enough to correlate the presently available dense gas 
with the star formation tracers. Note that the variation of the SFR 
cannot be too substantial: the time-integral of the current SFRs in 
the 100-pc ring and the 230-pc integrated CMZ over a Hubble time 
gives a total stellar mass that is roughly consistent with the total 
stellar mass enclosed at these radii (cf . Launhardt, Zylka & Mezger 
2002). 5 

If the star formation in the CMZ is episodic, it needs to be 
established which physics could be driving the variability. Could 
the low SFR simply be a stochastic fluctuation? At 0.015 Mq yr _1 
(Longmore et al. 2013a), the 100-pc ring produces about 10 5 Mq 
per dynamical time tdyn = 2n/Q, which corresponds to 2-3 young 
massive clusters (YMCs) and is thus consistent with the presence of 
the Arches and Quintuplet clusters. 6 If the SFR were consistently 
an order of magnitude higher (as predicted by density-dependent 
star formation relations), then ~ 25 such YMCs would be ex- 
pected. This would imply that the present cluster population in the 
CMZ is a <; 4.5a deviation. It is thus highly unlikely that the ob- 
served SFR is due to simple Poisson noise. 

Episodic star formation should be a common process in the 
centres of barred galaxies - in the inner Lindblad resonance rings, 
gravitational instabilities can drive the fragmentation of the nuclear 
ring and eventually induce a starburst, but this only takes place 
above a certain threshold density (Elmegreen 1994). If the system 
is steady-state on a global scale, this threshold density is given by 
Wburst = 0.6« 2 /G/imH and hence ?ib urst = {0.3,0.5,2.6} X 
10 4 cm -3 for the three regions of the CMZ listed in Table 1. This 
is comparable to the current mean density of no = 2 x 10 4 cm -3 , 
which may indicate that (part of) the CMZ is currently evolving to- 
wards a starburst. However, this critical density increases if the gas 
accretion rate along the bar is substantial. Additionally, orbital cur- 
vature causes the density waves in the central regions of galaxies to 
grow at an increasing rate towards smaller galactocentric radii, even 
in Toomre-stable discs or cases where the self-gravity of the gas 
does not set its global geometry (Montenegro, Yuan & Elmegreen 
1999). It is therefore easy to picture a system in which fresh gas is 
transported along the bar onto the 100-pc ring, where it builds up 
until the critical density for instability is reached - possibly at dif- 
ferent times throughout the ring. There is a notable population of 
24/im sources at / Ss 359.5° , i.e. beyond the position of Sgr C (see 
Figure 1 and Yusef-Zadeh et al. 2009), which may be the remnant 
of a recent, localised starburst. Variations would typically occur on 
the dynamical timescale of the system, which is ~ 5 Myr for the 
100-pc ring. 

Alternatively, feedback from YMCs might induce a fluctuat- 
ing SFR. As will be discussed in §3.1.3, the majority of the current 
star formation in the 100-pc ring takes place in the Sgr B2 complex. 
In the past, the birth environment of the 24/im sources (possibly 



3.1.2 Episodic star formation 

While it is tempting to assume that the CMZ is a steady-state sys- 
tem, the orbital and free-fall time-scales in the CMZ are so short 



5 Obviously, the bulge has not been in place for a Hubble time, but this 
comparison does put the present SFR in an appropriate perspective. 

6 We exclude the nuclear cluster of the Milky Way, of which the young stel- 
lar component is thought to have a very different origin (see e.g. Antonini 
2013 for a discussion). 
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also near Sgr B2) may have been the dominant site of star forma- 
tion. Such asymmetry implies that the feedback from YMCs origi- 
nates from discrete locations. This scenario is problematic, because 
the feedback energy will escape through the path of least resistance 
and hence YMCs on one side of the ring cannot support gas on the 
far side (also see §3.2.6). While a local starburst such as Sgr B2, 
or the one responsible for the formation of the 24pm sources, may 
not impact the entire 100-pc ring, it could blow out the sector con- 
taining the starburst. By the time such a burst is ~ 10 Myr old 
its Hll regions will have faded, and much of the gas could be in 
atomic form, being dispersed to higher latitudes by the combined 
acceleration of Hll regions, winds, and supernovae. 

The presence of a mere 2-3 clusters would necessarily place 
the current state of the ring at a minimum SFR, because the scarcity 
of feedback support suggests it should collapse to form stars over 
the next few free-fall times. During such a burst of star forma- 
tion, it could produce ~ 20 YMCs (which is the number of 
dense clumps observed by Molinari et al. 2011) of Arches-like 
masses, assuming that the progenitor clumps are similar to the 
cloud GO. 253+0.016 (also known as 'The Brick', see Longmore 
et al. 2012, with M ~ 10 5 M , R ~ 3 pc) and form stars in 
gravitationally bound 7 clusters at a 10-30 per cent efficiency. The 
timescale for the SFR fluctuations would be set by the half-time 
of the YMC feedback. Their bolometric luminosity decreases by 
a factor of two in roughly 8 Myr, implying that if the low SFR in 
the CMZ is explained by feedback-induced episodicity, a natural 
timescale would be ~ 10 Myr. 

A 5-10 Myr timescale for episodic star formation in the 100- 
pc ring places a strong limit on the observability of any produced 
YMCs - if their apparent present-day absence is not caused by 
their possible fading below the detection limit, their disruption 
timescale would have to be at most 10 Myr. The dynamical fric- 
tion timescale of a 10 4 Mq cluster orbiting in the 100-pc ring is 
tdf ~ 3 Gyr, 8 implying that dynamical friction does not affect 
YMCs on such a time-scale, and YMCs will therefore keep inter- 
acting with the ring regularly after their formation. The disruption 
of a 10 4 Mq cluster in a gas-rich, high-density environment is dom- 
inated by tidal shocks due to CMC passages (Elmegreen & Hunter 
2010; Kruijssen et al. 2011, 2012). The time-scale for disruption 
is then tf is oc Sq^ c p~^pymc (Gieles et al. 2006), with E G mc 
the GMC surface density, p mo i the spatially averaged molecular 
gas density, and pymc the YMC density. Taking the Brick as a 
template GMC, we have Egmc ~ 5 x 10 3 M Q pc~ 2 , whereas 
p m oi = £moi/2fr = 150 Mq pc -3 for the 100-pc ring (cf. Ta- 
ble 1). Substituting these numbers gives t s ^ B = 5.7 Myr for a 
10 4 Mq cluster with a half-mass radius of 0.5 pc (cf. Portegies 
Zwart, McMillan & Gieles 2010), which is much shorter than the 
*di^ al ~ 80 Myr due to the Galactic tidal field alone (Portegies 
Zwart et al. 2001). 9 The disruption of YMCs by tidal shocks is 



7 Using the cluster formation model of Kruijssen (2012) we find that at the 
high gas density of the 100-pc ring about ~ 50% of the stars are expected 
to form in bound stellar clusters. This idea is supported by the existence of 
the dispersed population of 24^tm sources (see Figure 1), and the roughly 
equal numbers of massive stars observed in clusters and in the field of the 
CMZ (Mauerhan et al. 2010). 

8 Using equation (30) of Antonini (2013), with 7 = 1.8, ro = r; n = 
80 pc, po = 145 Mq pc~ 3 , and m cl = 10 4 Mq. 

9 Note that Portegies Zwart et al. (2001) predict the presence of 50 Arches- 
like clusters in the CMZ based on a 80-Myr lifetime in the Galactic tidal 
field. When including tidal shocking due to GMCs, only three or four 
Arches-like clusters are expected to exist in the CMZ. Within the statisti- 



thus capable of 'hiding' the evidence left by previous starbursts. At 
2 and 4 Myr, the ages of the Arches and Quintuplet clusters are 
consistent with this disruption process. It therefore remains possi- 
ble that the star formation in the CMZ is episodic, which would 
then most likely be driven by local instabilities. 

Finally then, there might be a considerable time-variation of 
the gas inflow from large galactocentric radii onto the CMZ. Kim, 
Seo & Kim (2012) present numerical simulations of the gas flow 
in the central regions of barred galaxies, and measure the gas flux 
through a sphere with radius R = 40 pc, i.e. within their simulated 
equivalents of the 100-pc ring. They find that the variation of the 
gas flux is typically less than an order of magnitude in models with 
pronounced rings, because the inflowing gas is trapped in the nu- 
clear ring before gradually falling to the centre. For those models 
that do not develop such features, and hence have an uninhibited 
gas flow to the sphere where the flux is measured, the fluctuations 
sometimes reach two orders of magnitude. This can be taken as a 
rough indication of how the gas flow onto the 100-pc ring of the 
Milky Way may vary. The timescale for the variations in the gas 
flow corresponds to the dynamical timescale at the end of the bar, 
i.e. t dy n ~ 100 Myr, implying a relatively wide window during 
which the system can be observed at a low SFR. While this sce- 
nario therefore does not require us to put any stringent limits on 
YMC lifetimes, the preferred period for episodicity is shorter. The 
presence of two YMCs and the 24pm sources suggests substantial 
recent star formation activity, and considering the high densities of 
the several Brick-like clouds in the CMZ it not clear how long the 
current dearth of star formation will continue. 



3.1.3 The geometry of the CMZ and tidal shocks 

In highly dynamical environments like the CMZ or barred galaxies, 
it is conceivable that GMCs are disrupted by transient tidal pertur- 
bations ('tidal shocks', see e.g. Spitzer 1958; Kundic & Ostriker 
1995) before having been able to form stars. This idea was first in- 
troduced by Tubbs (1982), who suggested that perturbations would 
limit the star formation time-scale and hence decrease the SFR. 
Note that this effect differs from disruption by a steady tidal field 
(see e.g. Kenney et al. 1992 for extragalactic examples), to which 
we turn in §3.2. 1. The vast majority of the star formation in the 100- 
pc ring of the CMZ takes place in the complex Sgr B2 (Longmore 
et al. 2013a). It has been proposed by Molinari et al. (201 1) that the 
100-pc ring follows the X2 orbit - a family of elliptical orbits with 
semi-major axis perpendicular to the Galactic bar, which precesses 
at the same rate as the pattern speed of the bar, resulting in stable, 
non-intersecting trajectories. The X2 orbits are situated within the 
xi orbits, which are elongated along the major axis of the bar (see 
Athanassoula 1992). In this scenario, Sgr B2 and Sgr C lie at the 
points where the x\ and X2 orbits touch, and the accumulation of 
gas occurs due to two density waves with pattern speeds different 
than the flow velocity of the gas orbiting on the 100-pc ring. Such 
a configuration has been observed in other galaxies (e.g. Pan et al. 
2013). The ring upstream of Sgr B2 is fragmented into clouds that 
have properties suggesting that they could be massive protoclus- 
ters (Longmore et al. 2013b). A schematic representation of this 
configuration is shown in Figure 5. If Sgr B2 is indeed a stand- 
ing density wave and the impending encounter of these clouds with 
Sgr B2 is sufficiently energetic (cf. Sato et al. 2000), it is possible 

cal limits, this is consistent with the presence of the Arches and Quintuplet 
clusters. 
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Figure 5. Schematic representation of the picture discussed in §3.1.3 (also 
see Bally et al. 2010; Molinari et al. 201 1; Longmore et al. 2013b), as seen 
from Earth (top) and from above (bottom). Note that the position along the 
line of sight of the Arches and Quintuplet clusters is unknown, and therefore 
they are not included in the top-down view. As discussed in §2.1 and foot- 
note 2, it is possible that the ring extends further, in which case the x% orbits 
may connect to the ring under a different angle, and at a different galacto- 
centric radius. The discussion of §3.1.3 assumes the geometry depicted here 
to put upper limits on cloud-cloud interactions (see text). 



that star formation is briefly induced due to the tidal compression 
when they enter Sgr B2, but is subsequently halted when the clouds 
exit the Sgr B2 region and rapidly expand due to having been ti dally 
heated. We reiterate that the scenario of Figure 5 is not undisputed 
(see footnote 2). However, it does lead to the most extreme cloud- 
cloud encounters that could take place in the CMZ, and therefore 
we consider it as an upper limit. 

It is straightforward to quantify the disruptive effect of a tidal 
perturbation as the ratio of the energy gain AE to the total energy 
of the cloud E, under the condition that the duration of the per- 
turbation is shorter than the dynamical time of the perturbed cloud 
(the 'impulse approximation', see Spitzer 1987) - otherwise the in- 
jected energy is simply dissipated. In order to compute the relative 
energy gain, we consider a head-on collision and approximate the 
cloud and the perturber with Plummer (1911) potentials. If we in- 
clude the correction factors for the extended nature of the perturber 
(Gieles et al. 2006) and the second-order energy gain (Kmijssen 
et al. 201 1), the total relative energy gain becomes: 

AE 



1.1M^R-^\1 V 2 - 



(9) 



where Mr = M/10 7 M© is the perturber mass, Rh,i = -Rh/10 pc 
is its half-mass radius, ms = m/10 5 Mq is the cloud mass, 
rh,o = f*h/l P c is its half-mass radius, and V2 = V/100 km s _1 
is the relative velocity between both objects. The cloud is unbound 
if \ AE/E\ ^ 1. To describe Sgr B2 we adopt M 7 ~ 0.6 (Bally 
et al. 1988; Goldsmith et al. 1990) and R hA ~ 1, for the clouds 
approaching Sgr B2 we use the properties of the Brick, with m 5 ~ 
1.3 and rh,o ~ 2.1, and the relative velocity is taken to be V2 ~ 1 



(slightly higher than the line-of-sight streaming velocity in the 100- 
pc ring). For these numbers, equation (9) gives \AE/E\ — 2.8, 
suggesting that a Brick-like cloud could be unbound when pass- 
ing through the gravitational potential chosen to represent Sgr B2. 
To verify the validity of the impulse approximation, we note that 
the duration of the perturbation is At = 2R h /V = 0.1-0.2 Myr, 
whereas the dynamical time is tdya = (Gph) -1 ^ 2 ~ 0.4 Myr, and 
hence At < td yn - 

The above approach does not account for the detailed structure 
of Sgr B2 and the passing clouds, nor does it cover the collisional 
hydrodynamics of the clouds (e.g. Habe & Ohta 1992), their dissi- 
pative nature, or any possible, substantial deviations from spherical 
symmetry. However, it does make the interesting point that the en- 
ergy gain during an encounter of a Brick-like cloud with the Sgr B2 
region is comparable to the binding energy of the cloud. The clouds 
in the 100-pc ring are regularly spaced at ~ 25 pc intervals, and 
their evolutionary stage seems to be getting progressively older to- 
wards Sgr B2, with the clouds closest to Sgr B2 possibly having ini- 
tiated star formation, which is possibly caused by the triggering of 
cloud condensation at the orbital pericentre of the ring near Sgr A* 
(Longmore et al. 2013b). At a circular velocity of 80-140 km s _1 
(see §2.1), the Brick will enter Sgr B2 in 0.2-0.3 Myr and will 
leave it At = 0.1-0.2 Myr later. If the cloud is unbound by the 
interaction, star formation may only proceed for 0.3-0.5 Myr, or 
~ 2 global free-fall times. For a star formation efficiency (SFE) 
of 1-2% per free-fall time (Krumholz & Tan 2007), the clouds in 
the 100-pc ring will reach total SFEs of ~ 3%, which for the typi- 
cal cloud mass of 10 Mq implies that each cloud would typically 
produce 3000 Mq of stars. This is comparable to the total mass in- 
ferred from the 24/im point sources, and hence the above scenario 
may have applied to their formation too. Alternatively, the hydro- 
dynamic perturbation of the clouds passing through a density wave 
may also accelerate their collapse, in which case Sgr B2 and Sgr C 
represent the instigation points of star formation rather than the loci 
where it is terminated. 

We conclude that it is possible that the unique geometry of the 
CMZ plays a role in controlling star formation in at least part of 
the region. While this model may work well for star formation in 
the 100-pc ring, it remains to be seen to what extent similar effects 
could inhibit star formation elsewhere in the CMZ. The above esti- 
mate also ignores the possible contraction of the clouds and corre- 
spondingly shorter free-fall times during the passage of Sgr B2. We 
return to this scenario in §4.2, where possible tests are discussed. 

3.2 Locally inhibited star formation 

Here we discuss the possible mechanisms on spatial scales smaller 
than the disc scale height (AR < h) that may limit star formation 
in the CMZ. Volumetric star formation relations generally rely on 
the free-fall time ta oc n -1 ^ 2 to predict Esfr, in that approxi- 
mately 1% of the gas mass is converted into stars per free-fall time 
(Krumholz & Tan 2007). The relative universality of this number 
suggests that broadly speaking, there is always a similar fraction 
of the gas mass that participates in star formation. In star forma- 
tion theories, this is explained by the idea that the dispersion of the 
density PDF depends on the Mach number in a similar way to the 
critical density for star formation - the PDF broadens as the den- 
sity threshold increases (e.g. Krumholz & McKee 2005; Padoan & 
Nordlund 2011). Considering the broad range of gas densities ob- 
served in disc and starburst galaxies, this is a crucial ingredient to 
allow the SFE per free-fall time to be roughly constant. 

The above picture has been challenged by the recent obser- 
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vation by Lada, Lombardi & Alves (2010) that there is a possibly 
universal, critical volume density riLada ~ 10 4 cm -3 for convert- 
ing gas into stars. Such a transition should be expected, as there 
obviously exists some extreme density above which all gas ends 
up in stars (modulo the mass lost by protostellar outflows), but it 
is not clear why such a transition density would be universal. As 
discussed in §2, the fraction of the gas mass in the CMZ that is 
used to form stars is a factor of /sfr,u = 0.005 lower than the 
fraction of the gas with densities above the Lada threshold, in or- 
der for the SFR to be consistent with density-dependent (N > 1) 
star formation relations. We showed in §2 that if it exists, a typ- 
ical volume density threshold for star formation thus has to be 
nth ^ 10 7 cm~ 3 for various parametrizations of the density PDF, 
i.e. n t h = 10 7 -4 x 10 s cm~ 3 for a power-law approximation and 
n t h ~ {0.5, 2} x 10 9 cm~ 3 for a log-normal when including and 
excluding the effect of the magnetic field, respectively. 

In the following, we verify which physical mechanisms are 
consistent with the inhibition of star formation below such densi- 
ties. Because the turbulent pressure in the CMZ is remarkably high 
(e.g. Bally et al. 1988), with P tur b = no 2 ~ 10~ 6 erg cm" 3 , 
an important constraint is that potential star formation inhibitors 
should be able to compete with the turbulence. Therefore, we often 
use the turbulent pressure as a reference point to calculate the gas 
volume densities below which the SFR may be suppressed by each 
mechanism. The mechanism responsible for the low observed SFR 
needs to be effective up to a critical density of nth > 10 7 cm~ 3 
(see §2). 



3.2.1 Galactic tides 

A first condition for initiating star formation is that the progenitor 
clouds are not tidally disrupted. The gas needs to be at a density 
higher than the tidal density, i.e. the density required for a spheri- 
cal density enhancement to remain bound in a galactic tidal field. 
While this is not a guarantee that a region will eventually form 
stars, which requires time and space for it to become strongly self- 
gravitating, it does represent a key requirement for star formation 
to proceed. The tidal density is written as 



ntidal = 



3^4.pot Q 
47runiHG' 



(10) 



which only depends on the angular velocity fi. This expression 
assumes that the cloud orbit is circular. The constant ^pot de- 
pends on the shape of the galactic gravitational potential, and is 
^4pot = {0, 2, 3} for solid-body rotation, a flat rotation curve, 
and a point source (i.e. Keplerian) potential, respectively. Adopt- 
ing the Launhardt, Zylka & Mezger (2002) potential for the CMZ, 
we find A pot = {2.0, 1.9, 0.7} for the 230 pc-integrated CMZ, 
the 1.3° cloud, and the 100-pc ring, respectively. Galactic tides 
inhibit star formation in regions where ntidal > n t h, with n t h 
some unknown threshold for star formation. For the three regions 
of the CMZ that are listed in Table 1, equation (10) yields ntidal = 
{1.3, 1.8, 2.1} x 10 3 cm~ 3 . This is several orders of magnitude 
lower than the n t h ^ 10 7 cm~ 3 required to explain the SFR in 
the CMZ, but it is only an order of magnitude lower than the mean 
gas density in the CMZ. We did not account for eccentric orbits in 
this argument, because n t idai ^ n th requires such a high angular 
velocity (fi C rit ^ 10 2 Myr -1 ) that even at fixed circular velocity 
an eccentricity of e ^ 0.99 would be necessary for tides to be the 
limiting factor. We conclude that star formation in the CMZ is not 
inhibited by tides. 



3.2.2 Turbulence 

Turbulence plays a key role in the recent star formation models 
of Krumholz & McKee (2005) and Padoan & Nordlund (2011). 
While the former argue that turbulent pressure support sets the crit- 
ical density for star formation on the sonic scale, the latter work 
takes the point that turbulence is only responsible for driving lo- 
cal gravitational instabilities, and that the critical density for star 
formation is set by the thickness of the post-shock layers in the su- 
personic ISM. These differences aside, both models do predict a 
critical overdensity x = n/no for star formation 10 that scales with 
the GMC virial parameter a v ; r and Mach number as 



^turb = A x Ct v irM ■ 



(11) 



Both models also have remarkably similar proportionality con- 
stants A x ~ 1 to within a factor of 1.5. 

Adopting a virial parameter of a v ; r = 1.5 (cf. Krumholz & 
McKee 2005; Padoan & Nordlund 201 1), a Mach number of M = 
70 as in §2, and a mean density of n = 2 x 10 4 cm~ 3 , we see that 
the critical density for star formation is n tur b = :r tur bno ~ 1.5 x 
10 s cm -3 . This approaches the low end of the required threshold 
densities n t h ~ {0.5, 2} x 10 cm (see the dotted and dashed 
lines in Figure 3), although depending on the shape of the density 
PDF we indicated a strict lower limit of n t h ^ 10 7 cm~ 3 in §2. If 
we apply equation (11) to the Brick (Longmore et al. 2012), with 
Qvir ~ 1, M ~ 55, and no = 7.3 x 10 4 cm~ 3 , we obtain n tur b ~ 
2.2 x 10 8 cm~ 3 . By contrast, the typical properties of GMCs in 
the solar neighbourhood are a v ir ~ 1.5, M ~ 10, and no = 
10 2 cm -3 , which gives n turb ~ 1.5 x 10 4 cm~ 3 . Interestingly, 
this equals the Lada, Lombardi & Alves (2010) threshold for star 
formation to within the uncertainties of this calculation. 

The above numbers are very suggestive, in that the 'normal' 
overdensity thresholds for star formation on the one hand predict 
the Lada threshold for star formation in the solar neighbourhood, 
and on the other also predict a density threshold in the CMZ of 
n tur b ~ 2 x 10 8 cm~ 3 , which gives the best agreement with 
the observed SFR so far. The difference with the required nth = 
5 x 10 8 cm for a weakly magnetically confined density PDF is 
only a factor of three, 11 and may be offset by uncertainties in the 
density PDF, in the above numbers, and/or in the observed SFR. 
While turbulence is thus capable of increasing the density thresh- 
old to the required, extreme densities, we should note that this ar- 
gument is incomplete, as we have not established what is driving 
the turbulence. We return to this point in §5, and will now briefly 
discuss a particularly interesting uncertainty. 

3.2.3 A bottom-heavy initial mass function 

Observational measures of the SFR are strongly biased to the emis- 
sion from massive stars (m ib 8 M ). The low observed SFR in the 
CMZ may therefore be spurious due to an overproduction of low- 
mass stars with respect to the observed massive stars, whether the 
SFR is determined using massive YSOs (Yusef-Zadeh et al. 2009) 
or the ionising flux from massive stars (Murray & Rahman 2010; 
Lee, Murray & Rahman 2012). Neither technique is capable of re- 
liably detecting stars below 8 Mq . Recent observational studies of 
giant elliptical galaxies have found evidence for a bottom-heavy 

10 For the Krumholz & McKee (2005) model this assumes a typical size- 
linewidth relation for GMCs of a <x 5 . 

11 Note that the difference is still an order of magnitude when neglecting 
the (admittedly weak) effect of magnetic fields. 
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IMF An/ 'dm oc m " with a power-law slope of f3 — 3 at the 
low-mass end (m 1 Mq, see e.g. van Dokkum & Conroy 2010; 
Cappellari et al. 2012; Goudfrooij & Kruijssen 2013). 

It has been suggested that the characteristic mass scale of the 
core mass function (CMF) is set by the thermal Jeans mass roj 
(Elmegreen, Klessen & Wilson 2008) or the sonic mass m SO nic 
(Hopkins 2012). The Jeans mass is insensitive to the density and 
the radiation field if the Schmidt-Kennicutt relation is satisfied 
(Elmegreen, Klessen & Wilson 2008). Because in the CMZ it is not, 
the heating rate should be low compared to the pressure (relative to 
the Galactic disc) and hence the thermal Jeans mass is low too. This 
makes the CMZ a prime example of where the CMF could have 
a lower than normal peak mass. If we assume that the CMF and 
IMF are related, this suggests that the characteristic turnover of the 
IMF that is observed in the solar neighbourhood at m ~ 0.5 Mq 
(Kroupa 2001; Chabrier 2003) is environmentally dependent. Be- 
cause mj and m son i c decrease with the pressure, this results in 
an enhanced population of low-mass stars in the vigorously star- 
forming galaxies observed at high redshift (e.g. van Dokkum et al. 
2004; Daddi et al. 2007), which reach Mach numbers of M ~ 100 
(e.g. Swinbank et al. 2011) and may be the progenitors of current 
giant elliptical galaxies. Such Mach numbers strongly contrast with 
the M ~ 10 in the Milky Way disc. 

For core masses M ^ 1 Mq (i.e. stellar masses M Si 
0.5 Mq) and a Mach number of M ~ 100, Hopkins (2012) pre- 
dicts a steepened mass spectrum with a slope of ft ~ 3. Using pa- 
rameters that are appropriate for the Brick in the CMZ (Longmore 
et al. 2012) and adopting a SFE in protostellar cores of e — 0.5 
(Matzner & McKee 2000), we find em sonic ~ 0.01 M , whereas 
in the solar neighbourhood em son ic ~ 0.5 Mq. Similarly, at the 
approximate density threshold for star formation due to turbulence 
riturb ~ 2 x 10 8 cm~ 3 , the thermal Jeans mass in the CMZ 
(T = 65 K) is about emj ~ 0.06 Mq, whereas in the solar neigh- 
bourhood it is em j ~ 0.5 Mq. These low characteristic masses 
suggest that the CMZ might be the low-redshift equivalent to the 
progenitor environment of giant elliptical galaxies, and it is there- 
fore important to verify to what extent any unseen stellar mass at 
m < 0.5 Mq may increase the SFR inferred for a 'normal' IMF 
(Kroupa 2001; Chabrier 2003). 

We integrate the mass of a Kroupa (2001) IMF between 
m-mi-n = 0.08 Mq and m max = 100 Mq and compare it to the 
mass integral of a similar IMF, but with a power-law slope of ft = 3 
for masses m < 0.5 Mq. This increases the total mass by a factor 
of two, at the same number of massive stars. We conclude that while 
this is a non-negligible factor, it is (1) comparable to the uncertainty 
on the SFR in the CMZ and (2) in itself insufficient to explain the 
factor of ^ 10 suppression of the SFR in the CMZ. However, if 
we combine a doubled SFR with a density PDF that is narrowed 
by the high magnetic field strength in the CMZ (see §2.2), then 
the density threshold required by the observed gas depletion time- 
scale becomes nth ~ 2.5 x 10 8 cm~ 3 , which is comparable to the 
threshold predicted for normal turbulence (ri tur b ~ 2x 10 8 cm~ 3 , 
see §3.2.2). 

3.2.4 The atomic-molecular phase transition of hydrogen 

Dense molecular gas as traced by HCN is found to be corre- 
lated with star formation tracers (e.g. Gao & Solomon 2004; Wu 
et al. 2005). If this relation is causal in nature, molecular gas may 
be required to form stars (e.g. Schruba et al. 2011, although see 
Glover & Clark 2012 and Krumholz 2012 for an alternative view), 
implying that a star formation threshold could be related to the 



phase transition of Hi to H2 (Elmegreen & Parravano 1994; Schaye 
2004; Blitz & Rosolowsky 2004; Krumholz, McKee & Tumlinson 
2009b), which at solar metallicity occurs at E ~ 10 Mq pc~ 2 . This 
is the only theory that predicts a constant SFR per unit molecular 
mass (see §2). However, the density scale for the phase transition 
decreases with increasing metallicity (Krumholz, McKee & Tum- 
linson 2009a), and hence should be even lower than 10 M pc~ 2 
in the central bulge region (e.g. Brown et al. 2010). The bulk 
of the gas mass in the CMZ resides at higher surface densities 
and is indeed observed to be molecular (e.g. Morris & Serabyn 
1996). We therefore rule out the atomic-to-molecular transition at 
E 10 M pc~ 2 as the cause for a suppressed SFR in the CMZ, 
because the surface density is much higher in the central 1 kpc of 
the Milky Way. 

3.2.5 The Galactic magnetic field 

Another possible explanation would be that the high magnetic field 
strength in the CMZ (B <; 100^G, Crocker et al. 2010) inhibits 
star formation (e.g. Morris 1989). Using the condition that the mag- 
netic and turbulent pressure are balanced, this implies a critical den- 
sity of 

1 / B\ 2 

n mag = - , (12) 

2^oA tm H V J 

where [10 = 4n is the vacuum permeability constant. Using the 
values for each of the three CMZ regions from Table 1, this gives 
a critical density of n mag ~ 10 cm ~ 3 above which the magnetic 
pressure becomes less than the turbulent pressure. Therefore, the 
magnetic field cannot provide support against the turbulence in the 
CMZ. This result is unchanged when adopting the internal veloc- 
ity dispersions of the clouds in the CMZ rather than the large- 
scale velocity dispersion. However, it is important to note that the 
aforementioned value of the magnetic field strength applies to the 
low-density intercloud medium, and may be an order of magnitude 
higher in dense clouds (Morris, Uchida & Do 2006). This could 
increase the critical density to n mag ~ 10 4 cm~ 3 , which is com- 
parable to the mean gas density, but is still much lower than the 
density threshold required by density-dependent star formation re- 
lations. While it likely does not inhibit star formation, the presence 
of a 100/iG magnetic field is an important second-order effect. On 
the one hand, it is capable of narrowing the density PDF of the ISM 
in the CMZ somewhat, and slows down star formation accordingly 
(see §2.2 and §3.2.2). On the other hand, magnetic breaking leads 
to angular momentum loss during the condensation and contraction 
of cores, and hence accelerates star formation (Elmegreen 1987). 

3.2.6 Radiation pressure 

Considering the high angular velocity of the gas in the CMZ, 
the 100-pc ring needs to enclose some 10 9 M of mass. If we 
make the reasonable assumption that most of this mass is consti- 
tuted by stars, this implies a high stellar density of the CMZ of 
E* ~ 3 x 10 Mq pc~ within the gas disc scale height (see Ta- 
ble 1). It is therefore worth investigating whether stellar feedback 
is capable of inhibiting star formation in the CMZ. At first, feed- 
back is dominated by protostellar outflows, followed by radiative 
feedback, supernovae, and stellar winds. The relative importance 
of these mechanisms depends on the spatial scale and the environ- 
ment. It has been shown by Murray, Quataert & Thompson (2010) 
that in all but the lowest-density environments (e.g. GMCs in the 
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Galactic disc) radiation pressure is the dominant feedback mecha- 
nism for disrupting GMCs, whereas on scales <; 100 pc the energy 
deposition by supernovae becomes important. Crucially though, the 
total energy output from each of these mechanisms is comparable. 

We now test the hypothesis that stellar radiation inhibits star 
formation in the CMZ. As discussed in §3.1.2, this scenario has 
the problem that the CMZ hosts discrete star formation events, im- 
plying that the feedback on one side of the CMZ may not be able 
to affect gas on the opposite side. To test whether feedback is a vi- 
able explanation from an energy perspective, we again require pres- 
sure equilibrium between radiation pressure and turbulent pressure, 
which yields a critical stellar surface density of young stars S*, ra d 
above which radiation pressure is important: 



£*,rad — 



4-(l + </>trK0T 2 E) ' 



(13) 



where c is the speed of light, $ ~ 3 x 10 3 erg s" 1 g" 1 ~ 
1.5 x 10 3 Lq Mq -1 the light-to-mass ratio of a young, well- 
sampled stellar population (Thompson, Quataert & Murray 2005), 
and the term in parentheses indicates the optical depth 1 + <j) tr T 
with t ~ krT, ~ KoT 2 E, in which (f> t r = ftr/T ~ 0.2 
(Krumholz & Thompson 2012) is a constant that indicates the 
fraction of infrared radiation that is trapped at an optical depth 
r = 1, «r is the Rosseland mean dust opacity (cf. Thompson, 
Quataert & Murray 2005; Murray, Quataert & Thompson 2010), 
and Ko ~ 2.4 x 10" 4 cm 2 g" 1 K" 2 is a proportionality constant. 

We consider two cases for calculating the critical stellar sur- 
face density above which radiation pressure can compete with 
the turbulent pressure. Firstly, in the 100-pc ring we have no ~ 
2 x 10 4 cm" 3 , a ~ 35 km s" 1 and S = 3.0 x 10 3 M Q pc" 2 . 
Secondly, for the CMC 'the Brick' (Longmore et al. 2012) we 
adopt n ~ 7.3 x 10 4 cm" 3 , a ~ 16 km s" 1 and S = 5.3 x 
10 3 M Q pc" 2 . In both cases, we assume T = 65 K as in §2, which 

1. For the 100-pc ring, equation (13) 



1.6 x 10 5 M Q pc 2 . Combining this with the sur- 



yields an optical depth of r 
gives £* irad 

face area taken up by the gas in the ring (5 x 10 3 pc 2 ) and the 
lifetime of strongly radiating stars (~ 4 Myr, Murray, Quataert 
& Thompson 2010), we see that the 100-pc ring requires a SFR 
of ~ 200 Mq yr" 1 for radiation pressure to overcome the turbu- 
lence. This is four orders of magnitude higher than the measured 
~ 0.015 Mq yr" 1 . The Brick requires a similar surface density 
of E*, ra d ~ 1.1 x 10 5 M pc" 2 . Other than implying a criti- 
cal SFR of 0.7 M yr" 1 , it also means that a SFE of ~ 20 (not 
per cent!) is required to overcome the turbulent pressure, unless the 
newly formed stellar population has a five times smaller radius than 
its parent cloud, in which case a SFE of unity would imply similar 
turbulent and radiative pressures. The above numbers change by a 
relatively small amount when also including the combined flux of 
the more numerous, old stars of the Galactic bulge. Since radiative 
feedback is not capable of strongly affecting the gas on the cloud 
scale, it is also unable to stop the entire CMZ from forming stars 
(see 3.1.2). 



3.2.7 Cosmic rays 

Cosmic rays pressure could be important in the CMZ, either due to 
past star formation or black hole activity. For instance, it is possible 
that star formation in the CMZ is inhibited by the cosmic ray flux 
from supernovae (SNe). By equating the cosmic ray pressure to the 
turbulent pressure, and assuming the extreme case in which the cos- 
mic ray energy remains trapped in the CMZ, we obtain the critical 



volume density above which turbulence overcomes the cosmic ray 
pressure: 

"SN-EsNTsNTcr 

n-cr = TF ^ — • (14) 

In this expression, ??sn ~ 0.3 is the fraction of the SN energy 
that is converted to cosmic rays, Esn ~ 10 erg is the energy 
of a single SN, V is the volume, and Tsn is the SN rate, which 
for a Kroupa (2001) or Chabrier (2003) initial mass function (IMF) 
is given by r SN = 0.01 yr" 1 (SFR/M yr" 1 ), where the SFR 
is expressed in Mq yr" . The variable r cr oc n" 1 indicates the 
lifetime of cosmic rays to energy loss via collisions with protons in 
the gas, which is r cr ~ 10 s /no yr ~ 5 x 10 3 yr. 

Based on Table 1, we have SFRs of {3.3,0.06,1.5} x 
10" 2 M Q yr" 1 and volumes V ~ {3.3,1.4,0.05} x 10 6 pc 3 
for the 230 pc-integrated CMZ, the 1.3° cloud, and the 100-pc 
ring, respectively. Substitution into equation (14) then gives n cr = 
{0.1, 0.005, 3} cm" 3 (assuming a ~ 35 km s" 1 ), which for all 
three cases is much lower than the density threshold for star forma- 
tion required by the observed SFR. Given these numbers, it is very 
unlikely that cosmic rays from SNe affect the gas dynamics in the 
CMZ. For the cosmic ray pressure to compete with turbulence, the 
SFR would need to be four orders of magnitude higher than is be- 
ing observed. Even if all SN energy would be converted to cosmic 
rays, and their lifetime would be a thousand times longer, cosmic 
rays would still imply a critical density n cr < no. 

Alternatively, cosmic rays could originate from the activity of 
the central black hole of the Milky Way, which accretes at a rate 
of M <, 10" s M Q yr" 1 (Quataert & Gruzinov 2000; Baganoff 
et al. 2003). We assume that 0.5% of the accreted rest mass en- 
ergy is available to heat the gas (e.g. Di Matteo, Springel & Hern- 
quist 2005), and follow a similar argument as for the case of SN- 
powered cosmic rays above. This yields a critical volume density 
of tibh ~ 6 x 10" 3 cm" 3 for the accretion-generated cosmic 
ray pressure to compete with the turbulence in the volume of the 
100-pc ring. Cosmic rays are also unimportant compared to the 
thermal pressure anywhere other than the 100-pc ring, where the 
thermal and cosmic ray pressures are comparable to within a fac- 
tor of three. Hence, they may not be important kinematically, but 
they could be relevant for setting the temperature of the gas (Ao 
et al. 2013; Yusef-Zadeh et al. 2013). We conclude that feedback 
processes in general, and radiative, supernova, and black-hole feed- 
back in particular, are not consistently inhibiting star formation in 
the CMZ. 12 Of course, whether or not this also holds in other, ex- 
tragalactic cases depends on their recent star formation and black 
hole activity. 



4 IMPLICATIONS AND PREDICTIONS FOR FUTURE 
OBSERVATIONS OF THE CMZ 

We now turn to the implications and possible tests of the remaining 
plausible star formation inhibitors of §3, which are summarised as 
follows. On global scales, star formation could be episodic due to 
gas instabilities or variations in the gas inflow along the Galactic 
bar. Alternatively, the geometry of the CMZ could cause clouds to 

12 Although locally they may drive arches and bubbles - there will always 
be some (small) volume V such that the feedback pressure from the en- 
closed stars competes with the turbulent pressure. For the 1 00-pc ring (and 
hence Sgr B2 where most of the star formation takes place), cosmic rays 
overcome the turbulence for V Ss 8 pc 3 or R Ss 1 pc. 
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be heated by dynamical interactions. On local scales, the reduced 
SFR is consistent with an elevated critical density for star formation 
due to the high turbulent pressure in the CMZ. This solution would 
be aided by a weakly magnetically confined interstellar medium, 
which would narrow the density PDF, as well as by an underpro- 
duction of massive stars due to a possibly bottom-heavy IMF. 



4.1 Testing episodic star formation 

Evidence exists of episodic star formation events in the CMZ (So- 
fue & Handa 1984; Bland-Hawthorn & Cohen 2003; Yusef-Zadeh 
et al. 2009; Su, Slatyer & Finkbeiner 2010) and mechanisms have 
been proposed to explain how such episodicity can occur. As dis- 
cussed in §3.1.2, instabilities can drive the fragmentation of the nu- 
clear ring and eventually induce a starburst. la Gas in barred spi- 
ral galaxies like the Milky Way is funnelled from the disc through 
the bar to the galaxy centre (Sakamoto et al. 1999; Kormendy & 
Kennicutt 2004; Sheth et al. 2005). Because the conditions in the 
CMZ effect a higher threshold for star formation than in the disc, 
the gas needs time to accumulate before initiating star formation. 
While this by itself can already cause star formation to be episodic, 
it is also known from numerical simulations (Hopkins & Quataert 
2010) that the presence of a bar can cause substantial variations of 
the gas inflow towards a galaxy centre. The large variation of the 
central gas concentration of otherwise similar galaxies sketches a 
similar picture (Sakamoto et al. 1999). 

To constrain the possible episodicity of star formation in the 
CMZ, it will be necessary to map the structure of the gas flow along 
the Galactic bar, which is already possible using sub-mm and ra- 
dio surveys of the Galactic plane (e.g. Molinari et al. 2010; Walsh 
et al. 2011; Purcell et al. 2012). If the CMZ is presently near a 
low point of an episodic star formation cycle, then the gas needs to 
be accumulating and hence the inflow rate has to exceed the SFR. 
An improved 6D map of the CMZ itself would also help to un- 
derstand the nature of the possible, large-scale instabilities of the 
gas - the combination of line-of-sight velocities, proper motions, 
plane-of-the-sky positions, and X-ray light echo timing measure- 
ments should lead to a conclusive picture of gas inflow, accumula- 
tion, and consumption. Finally, the end result of the star formation 
process should be considered further. With infrared data and spec- 
tral modelling, it is essential to establish the recent (~ 100 Myr), 
spatially resolved star formation history of the CMZ. If evidence is 
found for a statistically significant variation of the SFR, the time- 
scale of such variations can be used to determine whether stellar 
feedback, gas instabilities, or a varying gas inflow rate are responsi- 
ble (although see §5.2-5.3 for additional constraints). Similar mea- 
surements could be made for nearby barred galaxies (see §5.4), 
paying particular attention to spectral features that signify a star- 
burst 50-100 Myr ago with no recent star formation, such as strong 
Balmer absorption lines and weak nebular emission lines (which is 
relatively straightforward using fiber spectra from the Sloan Digital 
Sky Survey, see Wild et al. 2007). 



3 Considering that the high Q gas (see §3.1.1) may inhibit the decoupling 
of gas clouds from the stellar background, we note that these are not likely 
to be gravitational instabilities, but are probably acoustic in nature (Mon- 
tenegro, Yuan & Elmegreen 1999). 



4.2 Testing interaction-limited cloud lifetimes 

The suppression of star formation by the unbinding of clouds in 
dynamical encounters has a simple observational implication: the 
gas volume density downstream of Sgr B2 should be lower than 
upstream, where the Brick-like (see Figures 1 and 5) clouds are 
currently residing. It is shown by Longmore et al. (2013b) that the 
gas density downstream of Sgr B2 indeed drops, but this could also 
be the result of star formation, and is only relevant to the problem at 
hand if the CMZ is a steady-state system. The option of interaction- 
limited cloud lifetimes (or their interaction-induced formation) thus 
requires a more detailed modelling in numerical simulations. 

A key test would be to see what happens to a Brick-like cloud 
when it enters the Sgr B2 region. A numerical simulation of the fate 
of the Brick in such a configuration should not be oversimplified - 
both the gas physics and the global dynamics have to be treated 
appropriately. If the Sgr B2 complex is indeed situated at the point 
where the x\ and X2 orbits coincide (similarly to the central re- 
gions in other galaxies, see e.g. Pan et al. 2013), the latter having a 
circular velocity of 80-140 km s _1 , the plausible interpretation is 
that Sgr B2 is a standing density wave - much like the spiral arms 
in the Milky Way disc. It should also be noted that the interaction 
between the co-moving gas flows on the xi and X2 orbits may very 
well influence the properties of the system. This means that the fu- 
ture evolution of the Brick is not accurately modelled by having it 
cross a fixed, gasless background potential chosen to reflect Sgr B2, 
not can it be covered by a simulation in which the Brick-like cloud 
is put on a collision course with a Sgr B2-like complex at the circu- 
lar velocity of the 100-pc ring. A more feasible approach would be 
to simulate the CMZ on a larger scale, including the gas flow along 
the bar. Such a simulation would be follow on existing work done 
on the gas flow in barred galaxies (e.g. Piner, Stone & Teuben 1995; 
Kim, Seo & Kim 2012), or on the type of large-scale, zoom-in sim- 
ulation of a disc galaxy that has the necessary spatial resolution to 
resolve the substructure of the gas on the X2 orbit (e.g. Hopkins & 
Quataert 2010). The key novelty of a new numerical investigation 
should be the production (or insertion) of Brick-like clouds on the 
X2 orbit and tracking their evolution in detail. We are aiming to 
address this aspect in a future work (Kruijssen, Dale, Longmore et 
al. in prep.). 

Observationally, it is not straightforward to infer the dynam- 
ical processes on the far side of Sgr B2, and the detailed study of 
nuclear rings is restricted to nearby (barred) galaxies. However, it 
is possible to map the properties of clouds like the Brick. ALMA 
data of the cloud kinematics and density structure will provide a 
key observational avenue to understand how such clouds may be 
affected by the high-density environment of the CMZ. 



4.3 Implications and tests of suppression by turbulence 

In §3.2.2, we show that a very promising explanation for the low 
SFR in the CMZ is an elevated density threshold for star forma- 
tion due to the high turbulent pressure (with M ~ 70), based 
on the star formation models of Krumholz & McKee (2005) and 
Padoan & Nordlund (201 1). This provides a theoretical justification 
for the conclusion made in §2 that the observed SFR is inconsistent 
with the universal threshold density from Lada et al. (2012). Even 
though the CMZ greatly extends the range of environmental con- 
ditions that can be probed in the Milky Way, a verification of the 
turbulence hypothesis requires our analysis to be extended to exter- 
nal galaxies. If the density threshold for star formation is as variable 
as would be expected on theoretical grounds, then there should also 
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be an environmental variation of which molecular line transitions 
correlate with star formation tracers. We intend to quantify this in 
a future work. 

Within our own Galaxy, ALMA observations of clouds in the 
CMZ are extremely well-suited for addressing the variation of the 
threshold density for star formation. This threshold needs to be 
higher in the CMZ than in the disc, and in the high-pressure en- 
vironment of the CMZ smaller spatial scales need to be probed in 
order to resolve the final stages of the cloud fragmentation and col- 
lapse towards cores and protostars. At the distance of the CMZ, 
only interferometers like ALMA and the EVLA are capable of 
reaching the necessary sensitivity and resolution in the wavelength 
range of interest. By combining star formation tracers with high 
density tracers, it will be possible to map the conversion of gas into 
stars in detail. A comparison with similar data from the solar neigh- 
bourhood should lead to a conclusive picture of the environmental 
variation of density thresholds for star formation. 

Note that for turbulence to be a conclusive explanation for the 
suppressed SFR in the CMZ, it is important to constrain the den- 
sity PDF in the CMZ - preferably well enough to verify the pos- 
sible influence of weak magnetic fields (see §2.2). It is possible to 
map the PDF by comparing the total flux above certain molecular 
line transitions. The one caveat is that these transitions will need to 
probe densities as high as the suggested threshold density of a few 
10 s cm -3 . It is also relevant to establish whether the IMF in the 
CMZ may be bottom-heavy due to the high Mach number and cor- 
respondingly low sonic mass. Current IMF determinations are only 
capable of reaching masses of m ^ 5 Mq (e.g. Bastian, Covey 
& Meyer 2010; HuBmann et al. 2012), but with the new genera- 
tion of large-scale facilities like the E-ELT it should be possible 
to probe the IMF at m < 1 M . A presently possible, but less 
direct method is to perform ALMA observations of the CMF in 
star-forming clouds in the CMZ. While this assumes a certain map- 
ping of the CMF to an IMF, the advantage is that star formation 
theories based on the sonic or thermal Jeans mass actually predict a 
CMF rather than an IMF, implying that a measurement of the CMF 
enables a more direct verification of our theoretical understanding 
of star formation in the CMZ. 



5 DISCUSSION 
5.1 Summary 

We find that the reduction of the star formation rate in the Cen- 
tral Molecular Zone of the Milky Way with respect to density- 
dependent (N > 1) star formation relations can be explained by 
the inhibition of star formation on both global and local scales. 

(i) Star formation could be episodic due to gas instabilities or 
variations in the gas inflow along the Galactic bar, with the present 
state of the CMZ corresponding to a near-minimum in the episod- 
icity. (§3.1.2) 

(ii) The geometry of the CMZ may limit cloud lifetimes by dy- 
namical heating, implying lower SFEs, but the phase-space distri- 
bution of gas in the CMZ is too uncertain to either confirm or rule 
out this scenario. (§3.1.3) 

(iii) Crucially, the high turbulent pressure in the CMZ increases 
the critical density threshold for star formation orders of magnitude 
beyond the Lada, Lombardi & Alves (2010) threshold, and hence 
substantially decreases the SFR with respect to such a threshold. 
(§3.2.2) To first order, this explains why the CMZ (Longmore et al. 
2013a) and its constituent clouds (Longmore et al. 2012; Kauff- 



mann, Pillai & Zhang 2013) are presently not efficiently forming 
stars. We discuss the possible mechanisms for driving the turbu- 
lence below. 

Interestingly, the density threshold implied by turbulent pressure 
for solar neighbourhood conditions is n ~ 10 4 cm~ 3 , which coin- 
cides with the Lada, Lombardi & Alves (2010) threshold (see §2.2 
and §3.2.2). This agreement theoretically supports the fact that a 
universal, fixed density threshold for star formation is ruled out by 
the observations. Second-order effects that could work in conjunc- 
tion with turbulence are a weak magnetic confinement of the ISM, 
which would narrow the density PDF (§2.2), as well as an under- 
production of massive stars due to a possibly bottom-heavy IMF 
(§3.2.3). 

A wide range of mechanisms is found to be unable to explain 
the dearth of star formation in the CMZ. 

(i) When including the stellar potential, the gas disc is found to 
be only marginally Toomre stable, and hence susceptible to gravita- 
tional collapse. While this rules out morphological quenching, the 
gas by itself is not strongly self-gravitating, which may slow the 
rate at which gas clouds can decouple from the stellar potential and 
become gravitationally unstable. (§3.1.1) 

(ii) The tidal density is lower than the mean density of the CMZ, 
indicating that star formation in not inhibited due to clouds being 
tidally stripped. (§3.2.1) 

(iii) The atomic-to-molecular phase transition of hydrogen oc- 
curs at too low a density to play a role in the CMZ, especially at 
high bulge metallicities. (§3.2.4) 

(iv) While the magnetic field is capable of affecting the density 
PDF of the ISM somewhat, the magnetic pressure in the CMZ is 
much lower than the turbulent pressure. (§3.2.5) 

(v) Radiation pressure should dominate the feedback energy at 
the surface density of the CMZ, but a SFE larger than 100% would 
be necessary to have radiation pressure compete with the turbu- 
lence. (§3.2.6) 

(vi) Cosmic rays, produced in supernovae or by the central black 
hole, are unable to overcome the turbulent pressure, even if the 
SFR was several orders of magnitude higher in the past. By con- 
trast, they do contribute to the thermal pressure in the 100- pc ring. 
(§3.2.7) 

The conclusion that feedback processes are not responsible for the 
low SFR in the CMZ also agrees with the discrete distribution of 
the main star formation events - if feedback did play a role, it would 
only be able to do so locally. 

It is interesting that the only global star formation relation that 
agrees with the low SFR in the CMZ is the Bigiel et al. (2008) 
relation (§2.1). Is this the 'fundamental' star formation relation? 
The mere existence of high-redshift starburst galaxies with Ssfr 
two orders of magnitude higher than the CMZ at the same gas sur- 
face density (Daddi et al. 2010b) suggests otherwise. As discussed 
throughout this paper, there are several reasons why the CMZ could 
be abnormal, and the extrapolation of the Bigiel et al. (2008) rela- 
tion to the CMZ therefore seems premature. The observational and 
numerical tests of the proposed mechanisms of §4 will be essential 
to draw further conclusions on this matter. 

In the remainder of this paper, we use the above findings to 
construct a self-consistent picture of star formation in the CMZ, as 
well as its implications for star formation in the centres of other 
galaxies. 
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5.2 What is driving the turbulence? 

While we find that (a combination of) more than one mechanism 
can be tuned to reach a satisfactory estimate of the SFR in the 
CMZ, it remains to be determined which mechanism is responsi- 
ble for inhibiting star formation. Out of the possible explanations 
listed in §5. 1 , the high turbulent pressure of the CMZ is perhaps the 
most fundamental. Contrary to the possible episodicity of the SFR, 
or perturbations acting on the cloud scale, it is the only mecha- 
nism that works instantaneously, at this moment in time. However, 
attributing the low SFR to turbulence alone is an incomplete ar- 
gument, because the turbulent energy dissipates on a vertical disc 
crossing time. Using the numbers from Table 1, the dissipation rate 
in the 230-pc integrated CMZ is roughly 5 x 1CT 21 erg cm~ 3 s _1 
(Mac Low & Klessen 2004). The dissipation time-scale is a mere 
£diss = 2h/a ~ 3 Myr, whereas for the 100-pc ring it is about 
0.5 Myr. Something is needed to maintain the turbulence. 

If we assume a steady-state CMZ, then the kinetic energy of 
the gas inflow along the bar is insufficient. The energy flux of the 
inflow is given by v 2 n fl/ii n fl, where v ln fi is the gas inflow speed, and 
£i n fl = Mgas/Minfl is the time-scale for accumulating the present 
gas mass. Under the assumption of a steady-state CMZ, the gas in- 
flow time-scale equals the time-scale for gas depletion t dcp \, which 
in the CMZ is about a Gyr (see Table 1). The energy dissipation 
rate is a 2 /t diss . Energy balance then gives 




(15) 



and substituting values representative of the CMZ yields a required 
Winfl for maintaining the turbulence of 600 km s _1 , which is much 
higher than the ~ 200 km s _1 that can be seen in the HOPS 
NH3(1, 1) data. If we relax the condition of a steady-state CMZ, 
a gas inflow time-scale of t- ln fi ~ 100 Myr (rather than the adopted 
^infl = 1 Gyr) would be needed to maintain the current turbulent 
pressure at the observed inflow speed. In order for such a scenario 
to be viable, the peak SFR would have to be more than ten times 
higher than is presently observed. This would put the CMZ tem- 
porarily on the star formation relations of equations (1) and (2). 

The gas inflow is only capable of driving the turbulence if star 
formation in the CMZ is episodic. In that case, the current gas de- 
pletion time-scale may not be an accurate estimate for the inflow 
time-scale. Specific combinations of gas inflow histories and star 
formation histories can explain the current degree of turbulence as 
the result of the global gas inflow. For instance, we might be seeing 
the product of a persistently high inflow rate, which builds up gas 
until finally a starburst is generated - with the CMZ presently at the 
stage between two bursts. 

Episodic starbursts would periodically drive superbubbles into 
the low-density ISM. As mentioned in §3.2.7, massive star feed- 
back can locally blow out parts of the 100-pc ring, producing the 
observed asymmetries in CO, NH3 and cold dust (see Figure 1). 
Radio and 7-ray observations of nested giant bubbles indeed pro- 
vide evidence that stellar feedback has a profound impact on lo- 
calized portions of the CMZ (Sofue & Handa 1984; Su, Slatyer & 
Finkbeiner 2010; Carretti et al. 2013). On large scales, the bubbles 
consist of swept up, low-density gas, whereas higher density mate- 
rial can also be driven to a bubble, but only on size scales AT? h. 

After cooling in the Galactic halo, these bubbles return to the 
disc as high-velocity clouds (e.g. Bregman 1980; Wakker & van 
Woerden 1997). Provided that the resulting energy flow is suffi- 
ciently high, this may drive the turbulence in the dense ISM. How- 
ever, we estimate in §3.2.7 that even at a star formation maximum, 



the energy density of feedback is probably smaller than the tur- 
bulent pressure. After conservation of energy (or a decrease due to 
dissipative losses) this implies that in-falling clouds are not likely to 
provide the necessary energy. While the return of ejected material 
therefore does not drive the turbulence, feedback itself is unlikely 
to do so for the same reason. Even if we imagine that feedback 
would be the source of the turbulence, then the presently low SFR 
implies that star formation would start after the turbulent energy 
has dissipated, i.e. in a few Myr from now. The resulting feedback 
would compress gas, which could trigger even more star forma- 
tion. Only if the gas were converted to a warm, diffuse phase, star 
formation could in principle be halted by feedback. However, the 
resulting blow-out would channel away the energy and a molecular 
disc would remain, which should still be forming stars. 

Quantifying the effects of feedback and other clas- 
sical turbulence drivers shows that neither are effec- 
tive. Following the compilation by Mac Low & Klessen 
(2004), we find that magnetorotational instabilities (yield- 
ing a heating rate of ~ 5 x 10~ 24 erg cm~ 3 s _1 ), 
gravitational instabilities (~ 10~ 24 erg cm -3 s _1 ), pro- 
tostellar outflows (~ 10~ 25 erg cm -3 s _1 ), ionizing 
radiation (~ 10 _2S erg cm -3 s _1 ), and supernovae 
(~ 2 x 10~ 23 erg cm" 3 s _1 ) are all incapable of compensating 
the dissipation rate. Only in the 100-pc ring, the heating rate due to 
gravitational instabilities (~ 10~ 21 erg cm -3 s -1 ) approaches the 
local turbulent energy dissipation rate (~ 10~ 20 erg cm~ 3 s _1 ). 

A more sustainable solution would be that the turbulence is 
driven by the (acoustic) gas instabilities themselves. As gas falls 
into the stellar-dominated potential of the CMZ, it grows unstable 
due to geometric convergence that compresses the gas to a higher 
density, even in the absence of self-gravity (Montenegro, Yuan & 
Elmegreen 1999), and the turbulent pressure increases accordingly. 
If the gas density has built up far enough that the compression 
leads to gravitational collapse, then the rate of energy dissipation 
increases and the gas eventually forms stars. 14 In this picture, the 
inner gas discs of galaxies are driven to a spiral catastrophe by 
acoustic instabilities, which drive the turbulence and produce irreg- 
ular gas structures 15 during an intermittent phase in which a minor 
role is played by the self-gravity of the gas and star formation. This 
naturally leads to bursty star formation, and once again it seems 
that even the explanation of the low SFR in the CMZ by turbu- 
lence requires some degree of episodicity. This scenario is consis- 
tent with the idea that the central rings of galaxies are prone to grav- 
itational instabilities when a certain density threshold is reached 
(see §3.1.2). 

Of course, there exists a global gas inflow rate above which the 
threshold density criterion for gravitational instability is always sat- 
isfied and star formation is no longer episodic. This may be one of 
the differences between the CMZ and the vigorously star-forming 
galaxies of Daddi et al. (2010b). The precise value of this critical in- 
flow rate varies with galaxy properties, but empirically it is straight- 
forward to distinguish both gas inflow regimes. If gas mass is build- 
ing up before undergoing a starburst, then the inflow rate must ex- 
ceed the SFR (even if the inflow rate is variable), whereas the SFR 

14 Note that the regions that reach the threshold density for gravitational 
collapse need not be randomly distributed throughout the CMZ, but may be 
corresponding to specific locations in the geometry, such as possibly Sgr B2 
(see §3.1.3 and Elmegreen 1994; Boker et al. 2008; Sandstrom et al. 2010). 

15 There is evidence for the presence of such dust lane structures in the cen- 
tral region of the Milky Way on larger (~ kpc) scales (see e.g. McClure- 
Griffiths et al. 2012 and references therein). 
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Figure 6. Schematic representation of a possible, unified picture of local and global processes that control star formation in the CMZ, as well as the conversion 
to observables and galactic star formation relations. See §5.3 for a more detailed discussion. 



of a starburst exceeds the inflow rate per definition. Conversely, if 
the inflow rate is so high that the gas mass does not have to build 
up before becoming gravitationally unstable and forming stars, the 
SFR is comparable to the inflow rate. This should also hold instan- 
taneously and would therefore be insensitive to variations of the 
inflow rate, as long as the inflow rate is high enough to persistently 
drive star formation. We therefore predict that the transition from 
episodic star formation to persistent, inflow-driven gravitational in- 
stabilities should occur around a mass inflow rate Minfl ~ Msfr. 
If the physics of the CMZ can be extrapolated to the high-redshift 
regime, then galaxies (or their centres) with M-mfi ^ Msfr should 
follow the Daddi et al. (2010b) sequence of galaxies with enhanced 
SFRs, whereas systems with Mi n fl > Msfr should statistically 
be found to have lower SFRs (making excursions to elevated SFRs 
and hence Mi n fl < Msfr during their periodic bursts of star for- 
mation). 16 

5.3 A self -consistent star formation cycle 

We have now discussed all the ingredients for a scenario in which 
the local and global explanations for the low SFR in the CMZ rep- 
resent different aspects of the same mechanism. This is crucial for a 
satisfactory explanation, because the different size scales influence 
each other. For instance, the overall gas supply is regulated on a 
global scale (AR > h), but a local condition must determine the 
gas consumption rate if turbulence sets the critical density for star 
formation. 

Based on the discussion in §5.1 and §5.2, we propose a multi- 
scale cycle that controls the star formation in the CMZ. It is sum- 
marised in Figure 6, which also illustrates the conversion of phys- 
ical quantities to the observables that part of the analysis in this 
paper is based on. The different stages of the cycle are as follows. 

(i) Gas flows towards the CMZ from larger galactocentric radii, 
which can be driven by secular evolution (e.g. dust lane or bar trans- 
port) or by external torques (e.g. due to minor mergers or gal axy 
interactions). 

(ii) Geometric convergence causes the inflowing gas to be com- 
pressed by acoustic instabilities, without the gas being gravitation- 
ally unstable. 

16 Because the detailed variation of galactic inflow rates over long time- 
scales is unknown, the present-day stellar mass (which gives a time- 
averaged SFR) can unfortunately not be used to constrain this argument 
further. For instance, the SFR of a system may presently be persistently high 
due to a high inflow rate, whereas in the past it could have been episodic at 
a low inflow rate. In that case, the time-averaged SFR is much lower than 
the current SFR, and combining the average SFR with the present inflow 
rate would incorrectly suggest current episodicity. 



(iii) The acoustic instabilities drive an increased turbulent ve- 
locity, which cascades down to smaller scales. 

(iv) The elevated turbulent pressure increases the local density 
threshold for star formation as nth oc noM 2 to n t h ~ 10 8 cm -3 . 

(v) Due to the high density threshold, the gas is not consumed to 
form stars and instead accumulates until the global density thresh- 
old for gravitational instabilities is reached. 

(vi) The gravitational instabilities drive the gas down to suffi- 
ciently high densities to overcome the local density threshold for 
star formation. 

(vii) At such high densities, the time-scale for turbulent energy 
dissipation is short, and the gas is rapidly turned into stars. This 
can occur at different times in different parts of the CMZ, as is 
exemplified by Sgr B2 undergoing a starburst while the other CMZ 
clouds are inactive. 

(viii) The cycle repeats itself, starting again from point (i). 

Whether or not the system is observed during a starburst or a 
star formation minimum depends on the relative time-scales of the 
stages in the above cycle. The inflow and accumulation of gas from 
stages (i)-(iv) take place on a galactic dynamical time-scale, which 
is 5 Myr on the spatial scale of the 100-pc ring, and ~ 100 Myr 
at the outer edge of the bar. As stated in §3.1.1, the time-scale 
for clouds to become gravitationally unstable during stage (v) is 
£ g rav ~ Qgas/ft (e.g. Jogee, Scoville & Kenney 2005), which in 
the 100-pc ring is t gra v ~ 1 Myr, but in the 230 pc-integrated CMZ 
it is tgrav ~ 20 Myr. This comparison shows that the 100-pc ring 
may be the unstable phase of the proposed cycle, and that stage (vii) 
may be reached at different times throughout the CMZ. Because the 
gas needs to reach high densities before becoming gravitationally 
unstable and forming stars, it is likely that much of the star forma- 
tion in the entire CMZ occurs in bound stellar clusters (see §3.1.2 
and Kruijssen 2012; Longmore et al. 2013b). 

The gas consumption time-scale during stage (vii) should be 
shorter than the time spanned by stages (i)-(v) on the scale of in- 
dividual clouds, because the free-fall time is only a few 10 5 yr or 
less at the high density of the clouds in the CMZ. However, the 
integrated starburst duration for the entire CMZ depends on mix- 
ing and hence should be comparable to its dynamical time-scale 
(i.e. up to 10 Myr), with the gas being resupplied on the dynami- 
cal time-scale of the bar (i.e. up to 100 Myr). Such relatively brief 
episodes of nuclear activity are thought to last up to several Myr, 
and are supported by HST observations of nearby galaxies (Martini 
et al. 2003). Statistically speaking, the CMZ should therefore be 
observed at a star formation minimum, although the simultaneous 
occurrence of mini-bursts of star formation in different parts of the 
CMZ could temporarily place it on or above galactic star formation 
relations. 

Verifying our proposed scenario using galactic scaling rela- 
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tions alone is not straightforward. One qualitative constraint fol- 
lows from Figure 2, which shows that on a certain length scale 
(AR = 230 pc) the CMZ fits the empirical Schmidt-Kennicutt 
relation of equation (1). It may be that on such a spatial scale, all 
of the stages (i)-(vii) are being sampled. However, even when av- 
eraged over 230 pc, the CMZ does not fit the Silk-Elmegreen re- 
lation of equation (2). Krumholz, Dekel & McKee (2012) recently 
proposed that a universal star formation relation is obtained by sub- 
stituting — > l/t s f in equation (2), where the star formation time- 
scale t s f is either set by the local (GMC) or global (galactic) dy- 
namical time-scale, whichever is the shortest. If we account for the 
long condensation time-scale of GMCs and hence increase t si by a 
factor of Qgas (as in the Toomre regime of Krumholz, Dekel & Mc- 
Kee 2012), we see that the 230 pc-integrated CMZ fits a modified 
Silk-Elmegreen relation Esfr = ^SES^/Qgas- The 1.3° cloud 
and the 100-pc ring do not, and in the framework of the cycle of 
Figure 6 this is because on such small scales the episodicity of star 
formation is important. 

The conversion of the proposed cycle to observables is well 
beyond the scope of this paper, but Figure 6 shows a qualitative il- 
lustration. The global gas surface density mainly samples stage (v), 
whereas the star formation rate density traces stage (vii) of the cy- 
cle. However, the relation to the actual volume densities of star- 
forming gas and young stars depends on several factors. Obser- 
vational limitations such as the spatial resolution, the sensitivity, 
and possible tracer biases lead to measurement uncertainties in the 
derivation of the relevant physical quantities. For instance, unre- 
solved stars and gas may not be occupying the same volume (e.g. 
due to different filling factors or porosities, Silk 1997), and the sam- 
pling of the gas (PDF) and the star formation may be incomplete 
(the degree of which depends on the adopted tracers). Finally, there 
may be systematic conversion biases (e.g. Xco) when translating 
measured quantities to physical ones. Before a consensus on these 
points exists, it is more feasible to verify the proposed scenario as 
suggested in §4. 1 and §4.3, i.e. focussing on the parts rather than the 
sum thereof. For example, stages (i)-(iii) may be probed in nearby 
galaxies, with global gas flow kinematics constraining stages (i) 
and (ii), and velocity dispersion measurements tracing stage (iii). 

5.4 The extragalactic regime 

Given that the CMZ can only be observed in its present state, the cy- 
cle of Figure 6 is sampled better by including other galaxies where 
the same process is taking place. The inner regions of galaxies of- 
ten have dust structures with the number of spiral arms increasing 
with radius, and power spectra that are consistent with them being 
a manifestation of acoustic turbulence (e.g. Elmegreen et al. 1998; 
Martini & Pogge 1999; Elmegreen, Elmegreen & Eberwein 2002). 
These structures host little to no star formation, which if present is 
often arranged in a circumnuclear ring (e.g. Barth et al. 1995; Jogee 
et al. 2002; Sandstrom et al. 2010; van der Laan et al. 2013), and 
their gas by itself is often not strongly self-gravitating (Sani et al. 
2012). In all respects, the central regions of these galaxies are rem- 
iniscent of the CMZ. Indeed, there are known examples of nearby 
galaxies of which the star formation activity in their central regions 
falls below galactic scaling relations (Hsieh et al. 2011; Nesvadba 
et al. 2011; Sani et al. 2012), even in cases which at first sight may 
be perceived as a starburst system (Kenney, Carlstrom & Young 
1993). By contrast, examples of galaxy centres with a normal or 
enhanced SFR exist too, which may indicate a different stage in 
the same process (e.g. the central starbursts in the sample of Jogee, 
Scoville & Kenney 2005, also see Sakamoto et al. 2007, 2011). A 



statistical census of the activity in such galaxies enables a quantita- 
tive estimate of the duration of the different stages in Figure 6. This 
is a key step, which has been made previously. For instance, Martini 
et al. (2003) consider a sample of 123 nearby, active and passive, 
barred and unbarred galaxies, and estimate that the quiescent stage 
takes up to ten times longer than the star formation episode, which 
is consistent with what we propose for the CMZ in §5.3. 

Following a similar argument, Davies et al. (2007) find that 
active galactic nuclei (AGN) typically become active some 50- 
100 Myr after the onset of prior (but ceased) starburst activity. 
This is consistent with the scenario of Figure 6 - a starburst and 
its corresponding feedback may either consume or eject the gas, 
preventing the black hole from being fed until supernova feedback 
has ended (after 40 Myr). After that time, the gas flow towards the 
black hole continues. Because in our scenario the gas is not (yet) 
self-gravitating, it can proceed to the black hole without forming 
stars. Due to the gas inflow and preceding starburst activity, this 
phase should be characterized by centrally concentrated surface 
brightness and gas mass profiles, which has indeed recently found 
by Hicks et al. (2013). We conclude that our results may well have 
implications for the massive black hole growth and AGN activity. 
Specifically, the phase of low star formation activity in the cycle of 
Figure 6 could provide a window for efficient black hole growth. 

One of the key questions is how the physics that regulate star 
formation in the CMZ and the central regions of nearby galaxies 
translate to other extragalactic environments. An interesting target 
would be giant elliptical galaxies that contain a sizable, but remark- 
ably quiescent gas reservoir (e.g. Crocker et al. 2012). In terms 
of the gas properties, the CMZ is remarkably similar to the con- 
ditions found in ULIRGs and high-redshift, star-forming galaxies 
(Kruijssen & Longmore 2013). However, these galaxies are under- 
going vigorous starbursts, whereas the SFR in the CMZ is much 
lower than expected. We suggest in §5.2 that rapidly star-forming 
galaxies may have such a high global gas inflow rate that their con- 
stituent gas can always become gravitationally unstable, rather than 
having to build up to a certain density and forming stars episod- 
ically. This hypothesis can be tested with a simple observational 
prediction. Galaxies or their centres with Mind Ss Msfr should 
follow the Daddi et al. (2010b) sequence of galaxies with enhanced 
SFRs, whereas systems with M- ln a > Msfr. should statistically 
be found to have lower SFRs (making excursions to elevated SFRs 
during their periodic bursts of star formation). While observation- 
ally challenging, high-resolution ALMA observations should en- 
able this prediction to be tested (see §5.2 for a more detailed discus- 
sion). Alternatively, the star formation in a subset of such galaxies 
may take place at sufficiently large galactocentric radii for acous- 
tic instabilities to be unimportant, 17 i.e. in the regime where the 
gravitational instability of the gas alone is driving star formation. 
Another obvious concern is that by specifically considering rapidly 
star-forming galaxies, we are already biasing ourselves to what may 
be an episodic peak of star formation activity - in which case star 
formation in these galaxies could actually be similar to that in the 
CMZ over longer time-scales. 

In summary, we find that the lower-than-expected SFR in the 
CMZ can be explained by a self-consistent cycle connecting the 
galaxy-scale gas flow, acoustic and gravitational instabilities, tur- 
bulence, and local star formation thresholds. The cycle is summa- 
rized in §5.3 and Figure 6. We conclude that a low SFR may well 
be a common property of the centres of Milky Way-like galaxies, 



The circumnuclear starburst of Arp 220 is a clear counterexample. 
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although examples of the opposite should exist too if our proposed 
explanation is valid. The question remains how the lessons learned 
about star formation in the central regions of galaxies can be extrap- 
olated to other galactic environments. As discussed in this section 
and §4, the combination of future numerical work with data from 
new observational facilities should allow this topic to be addressed 
in the necessary detail. 
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